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ABSTRACT 


Segregation Roasting process for oxidised Copper and 
Nickel Ores has been investigated. Both the theoretical and 
experimental aspects have been •xamined. The thermodyn.amic 
analysis of the Segregation system involves the solution of a 
system of simultaneous nonlinear algebraic equations. In the 
present study, two numerical techniques namely Brown’s method 
and Powell’s hybrid method that overcome the deficiencies in the 
classical methods of solution of these equations have been 
employed. It is observed that both the methocta lead to fast 
ultimate convergence, although Powell’s method can accomodate 
a larger span of initial estimates of the unknowns. Solution 
of the equations has been used to assess the feasibility of the 
process in the case of Cu and Ni Segregation, and to predict 
the effect of the process parameters like temperature, amounts 
of chloride, reductant and moisture. The analysis indicates 
that qualitatively. Segregation is favourable in the case of 
Copper while some serious constraints hinder Nickel Segregation. 
Coupled with this, a generalised computational procedure that 
utilises the principle of minimisation of Gibbs' free energy 
has been attempted. It provided the equilibrium composition 
of both the gaseous and condensed species. 

Laboratory scale experiments on Segregation Roasting has 
been carried out to atleast qualitatively predict the operating 
conditions required for Segregation of synthetic Oxidised Copper 
Ore. 
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CHAPTER 1 


INTRODUCTION 

Oxidised Copper Ores are commonly concentrated either 
by Flotation or by Leaching. However, flotation is not 
possible when the copper minerals are locked with the gahgue 
or when' the copper minerals are mainly silicates, and leaching 
is uneconomic when the presence of limestone or dolomite in 
the gangue causes excessive acid consumption. 

The segregation roasting process (which now onwards would 
be referred to as segregation process) provides' an alternate 
method for treating oxidised copper ores by extracting the 
copper from the ore minerals at red heat and reprecipitating 
it as a metal on finely divided carbon which is admixed with 
the Ore. After cooling, in the absence of air, or quenching, 
the copper is recovered in a high grade concentrate by • 
flotation. 

In practice, the process involves the addition of 
1-1. fine coal or coke and 0,3 to 1 percent sodium chloride 
to the ground Ore, and heating to temperature of the order of 
800°C. The basic steps considered to occur are ; 

1. The salt is hydrolysed by moisture in the presence of 
silica, Aluminosilicates and other clay minerals to 
form hydrogen chloride. 
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2. The Hydrogen chloride reacts with the Copper Oxide in the 
Ore producing H 2 O and a volatile chloride, generally 
agreed to be predominantly the copper (I) chloride trimer, 

3. The Cuprous chloride, is reduced by Hydrogen in the vicinity 
of the carbon particles to form the segregated metallic 
copper, at the same time regenerating Hydrogen chloride. 

Under optimum conditions mild reducing conditions prevail 
in the segregation chamber, and thus the copper concentrates 
in the stronger reducing environment of the carbonaceous 
reagents. • 

The process continues in a cyclic manner and amounts of 
chloride and Hydrogen considerably less than the stoichiometric 
are required. The segregation process also recovers Gold or 
Silver that may be in the Copper Ore - metals that are usually 
lost in the tailings during Acid or Ammonia Leaching. Ores of 
many metals that form volatile chlorides or Oxy-chlorides may 
respond to the segregation process. It has been suggested that 
Antimony, Bismuth, Cobalt, Gold, Lead, Nickel, Palladium, 

Silver and Tin may be treated in this manner. In recent years, 
the process has been extended to roasted concentrates and 
selective segregation of Copper from Copper-Nickel concentrates. 

Scope of Present Work ; 

In our country, there are deposits of oxidized copper and 
Nickel Ores which have hitherto not been exploited fully due to 
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the lack of a viable process for extraction, and in this 
background, the segregation process seems to be a potential one 
that necessitates a deeper theoretical and experimental study. 

As elucidated earlier, segregation process is typical of those 
systems involving many or complex equilibria, a thermodynamic 
evaluation of v^hich involves the numerical solution of a system 
of simultaneous nonlinear equations. These are generally solved 
by : 

1) Trial and Error approach or graphical method 

2) Classical methods of solution of nonlinear equations. 

The first approach becomes extremely tedious and time consuming 
when many sets of numerical values of the thermodynamic varia-r 
bles are required? the second approach not very infrequently 
runs aground when the variables vary over a wide range, or 
when the flat start values are far away, from the actual zeros 
of the oquations. 

In the present study, Algorithms that overcome these 
deficiencies are applied to assess tho thermodynamic feasibility 
of the segregation process for the extraction of Cu and Ni from 
low grade Oxide Ores. Further, the effect of a number of process 
variables on segregation is predicted. Although tho numerical 
technique is exemplified for the segregation process, it would be 
applicable for any heterogeneous system. Further, a generalised 
computational scheme that includes both the gaseous and condensed 
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species has been employed in the determination of equilibrium 
compositions by minimising the Gibb’s Free energy. Secondly, 
laboratory batch tests on synthetic Copper Ore for the 
extraction of Copper by the segregation process were designed to 
study the effect of process variables. 
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CHAPTER 2 

LITERATURE REVIEW 

The subject of Segregation Process has spawned 
an enormous literature under various headings ; Segre- 
gation chemistry^ Thermodynamic and Kinetic 
aspects^^^ Pilot and Commercial Operations^ , 

Effect of various Reagents on Segregation etc., an exhaustive 
review of which is given by Habashi^^^. In this review, only 
the Thermodynamic aspects of Segregation Process are consi- 
dered. 

Hitherto, a detailed study from the Thermodynamic stand- 
point has been made by Iwasaki^^^ and Yazawa^^^ I the former 
by employing Oxygen-chlorine potential diagram, and the latter 
by combining stoichiometric relations with equilibrium con- 
stants . 

Iwasaki^^^ constructed Equilibrium diagrams showing the 

stability regions of oxides and chlorides of metals, 
and superimposing them for comparison. 

An equilibrium diagram for Cu-H^-HCl and Fe-H 2 -HCl at 750°C 
and at constant partial pressures of v^rater vapor of 0*3 atm 
and 0.03 atm were constructed. It is found that the effect 
of water vapor concentration on the relative positions of the 
equilibrium lines is small for both copper and iron. There is 
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a marked difference in the stability fields of the 2 systems 
indicating that the Copper Oxides may be readily converted to 
CuCl, whereas the Iron Oxides remain as Oxides. The figure 
also indicates i 

1. that too strong a reducing atmosphere may reduce the 
copper oxides to metallic copper in place, and the result- 
tant metallic copper may not be chlorinated (reduction 

in situ) 

2. that a moderately reducing atmosphere helps the chloridi- 
sation of Cupric Oxide 

3. an excessive amount of Oxygen in the atmosphere would 
interfere with the process by preventing the conversion of 
both CuO and CuCl and also of CuCl to metallic Cu* 

A similar diagram involving the systems Ni-H2“HC1 and 
Fe-Iip-HCl at 950°C and at the water vapor concentrations of 
3'/. and 30/. by volume were constructed. The stability fields 
of Nickel and Iron chlorides overlap very closely and the 
partial pressure of HCl required to chloridize Nickel is a 
few orders of magnitude greater than that for Copper which is 
in good agreement with a laboratory observation that a rela- 
tively large quantity of Salt is required to effect Nickel 
Segregation and that CaCl 2 is more effective than NaCl. 
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A metal chloride formed must be volatilised and transv 
ported to a precipitation site, most commonly at the surface 
of a solid reductant where it is reduced. Brewer and Lofgren^^^ 
established that cuprous chloride exists as the trimer 
Cu^C^l^lg) v/hich dissociates into the monomer CuCl(g) at high 
temperatures 5 that cuprous chloride is almost entirely in the 
form of the trimer in the neighbourhood of 8CX3°C, and that no 
evidence was found for the presence of the dimer. 

(27) 

Kubaschewski and Evans^ ' plotted the vapor pressures of 

the chlorides of several metals of interest in segregation 

roasting against temperature. A high equilibrium vapor 

pressure may be indicative of a successful transport' of a metal 

(l9 20) 

chloride. According to Rey'^ '^Ferrous chloride has a high 
enough vapor pressure, but Iron cannot be segregated along with 
Copper since water vapor decomposes Iron chlorides to form 
Magnetite. The vapor pressure of AgCl in the relevant tempera- 
ture range being very low (10 " mm Hg), Segregation of Ag is 
unsuccessful. Apparently, Copper functions as a carrier of 
Gold and Silver. 

It is observed that each element has the relatively well 
defined regions of stability of its Oxides and chlorides and 
the relative positions of these regions change very little with 
temperature or water vapor concentration. This lead Iwasaki^ 
to construct equilibrium diagrams at 750*^0 of several metals of 
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interest in Segregation Roasting . From these diagrams, it may 
be speculated that Pb may respond to Segregation;. Iron and 
Manganese may not respond to Segregation as their Oxides are 
the stable species (at this temperature). 

According to Rey^^*^’^^), segregation atmosphere cleans the 
gold surface and improves its flotability. Cobalt appears to 
be more easily chloridisod than Nickel, but the partial pressure 
of hydrogen required to reduce it is higher than that required 
by Nickel. 

The accepted mechanism on the precipitation of metallic 
copper on carbon particles is the reduction of cuprous chloride 
vapor by in the vicinity of carbon particlesf^^^ 

CU3CI3 + 3/2 H2 = 3Cu + 3HC1 

Inability of carbon as a reducing agent for metal chloride may 
be inferred from the thermodynamic consideration of the reaction 
(Kellog^^^ 

2MCI2 + C = 2M + CCI4 

( 20 ) 

Rey ' j however, showed that in the presence of a small amount 
of water vapor, its reaction with carbon produces hydrogen, and 
the reduction of metal chlorides becomes possible. The presence 
of a solid reductant that provides nuclei for precipitation and 
growth of metallic copper from cuprous chloride vapor, and that 
releases H2 at a controlled rate,, appears essential for the 
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process. Thus less active the carbonaceons material, the 
more difficult it will be for Copper to precipitate. 

A thermodynamic Evaluation of Segregation process was 
( 2 ) 

performed by Yazawa^ * by assuming establishment of local 
equilibria in each step at each site. This was done by 
juxtaposing the thermodynamic relations arising out of the 
Mass action equations with the stoichiomatic relationships 
for which, feed pellets of Copper and Nickel Oxides with the 
typical compositions were segregation roasted in a Lab scale 
Experiment. The solution of the resulting set of nonlinear 
equations were presented in a plot of the partial pressures 
of the various gaseous species against Oxygen potential. From 
the resulting profiles, of Copper and Nickel Segregation 
systems, Yazawa concluded that the fundamental reaction cal- 
culations are promising for Copper whereas serious restrictions 
hinder Nickel Segregation. 
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CHAPTER 3 


STRUCTURE OF THE PROBLEM 


The Segregation Process involves a complex set of 
reactions, occurring at adjacent sites in a mixture consisting 
of the Ore, Salt and Carbon, For a thermodynamic analysis, 
chemical equilibrium is assumed in the system, and the equili- 
brium constant method is employed. 


In a typical Cu Segregation system, the predominant gaseous 
species at equilibrium are CO, CO 2 , H 2 » H 2 O, CuCl, Cu^Cl^, HCl 
and FeCl 2 . For the 3 nonmetals C , H, Cl, expressions for the 
conservation of masses may be written, 'nji^* indicating the 
number of moles of the chemical aspecies i. These expressions 
are : 


C balance 
H 2 balance 


= "c = "CD "OOg 

: "H2= "112 + "HjO + 1/2 "h:i 


(3.1) 

(3.2) 


Cl balance 


: '^Cl = ”hC 1 + ^“^FeCl2 + ^”Cu3Cl3 + ”CuCl 

(3.3) 


the partial pressures of the gaseous species i are defined by 

Hi 

pi 

being the total number of moles, and P, the total pressure. 
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Here , 


Ht = ’^HCl + n. 


CuCl '^Cu3Cl3 ”FeCl2 

( 3 . 4 ) 


Equations ( 3 . 5 ) to ( 3 . 11 ) constitute the equilibrium identified 
in the Copper Segregation system along with the expressions 
for the equilibrium constant at a particular temperature. 

CU20(S) + 2 HCl(g) = 2/3 003013(9) + H20(g) 


2/3 

^5 ‘^ 003012 * 


H 2 O* '^T 


l/ 3/^2 


HCl 


( 3 . 5 ) 


2/3003013(9) + H2(g) = 2 Cu(S) + 2 HCl(g) 
^6 *‘^H2 • 


3 3 


( 3 . 6 ) 


Fe203(S) + 4Hcl(g) + H2(g) = 2Fecl2(g) + 3H20(g) 

^34 

S ^®^^2*'^H2 o/’^HC1 • ’^H2 
Fe 304 (S) + 6HCl(g) + ^(g) =3FeCl2(g) + 4H20(g) 


y. 3 4 /y>^ n 

^8 “ ^FeCl^^H20''^HCl H 2 


( 3 . 7 ) 


( 3 . 8 ) 


1/2 Fe2Si04(S) + 2 HCl(g) = FeCl2(g) + 1/2 Si02(S) + H20(g) 


^9 "" ’^Feei2’‘^H20/'^Hei 


( 3 , 9 ) 
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CO(g) + H20(g) == 002(9) + H2(g) 


'^10 “ ''cD.2-"H2^"cX)-'’H20 
3CuCl(g) = Cu2Cl2(g) 

^11 ^ *^003013 *’^t/^CuC1 


(3.10) 


(3.U) 


If the value of si'e known (from experiments) 1 the 

following equations may be solved simultaneously, at a given 
CO 2 /CO ratio : 


Eq. (3.1), (3.2), (3.3) and (3.4) 
Eq. (3.5) or (3.6) 

Eq. (3.7) or (3.8) or (3.9) 

Eq. (3.10) or (3.11) 


The equilibrium (Pqq /p 


qq) ratios at which the solid species 
Cu and CuO are stable is delineated in 


Appendix 3.2* 


The solution of the system of simultaneous Nonlinear 
Algebraic equations outlined above yield the localised concen- 
trations of the various gaseous species. 


In a subsequent chapter, the feasibility of the process 
is assessed from the plot of the partial pressures of the 
gases generated against the Oxygen potential. 
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The stoichiometric values may be varied from the base 
value, and comparison of the resulting partial pressures of the 
various gaseous species with the base value aids in the predic~ 
tion of the effect of the process variables in the segregation 
system. 

Nickel Segregation Syst em 

Analogous treatment as in the case of Copper may be extend 
ed to Nickel for a thermodynamic evaluation. In the Nickel 
Segregation system, the predominant gaseous species constitute 
CO, CO 2 , H 2 > H 2 O, Hcl, NiCl 2 and FeCl^, 

The stoichiometric relations are s 


C balance 
H 2 balance 
CI 2 balance 


• ’^H2 ~ ”H2 

* '^Cl2 “ ’^HCl ’^FeCl2 ’^NiCl2 

"^T ’^CO “^CC^ ^ 

“^HCl ’^NiCl2 '^FeCl2 


(3.12) 

(3.13) 

(3.14) 

(3.15) 


The thermodynamic relations are i 

NiO(S) + 2HCl(g) = NiCl 2 (g) + H20(g) (3.16) 


% “ ^N iCl 


2* ^H20'^^HCi 
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NiCl2(g) + H2(g) = Ni(S) + 2HCl(g) 


(3.17) 


^6 “ ’^HCl/'^NiCl^* 

Fe304(S) + 6HCl(g) + H^Cg) = 3FeCl2(g) +4H20(g) 

^ _ 3 4 ; 6 

^7 ■■ ^FeCl2‘’^H20 /'^HC1*'^H2 


(3.18) 


1/2 Fe2SiO^(S) + 2HCl(g) = FeCl2(g) + 1/2 Si02(S) + H20(g) 

(3.19) 

^8 "" ’^FeCl2'”H2o/’^HCl 


CO(g) + H20(g) = C02(g) + H2(g) 


(3.20) 


S = "c02-''H2/''c0-"H20 


As in the case of copper, outlined in the previous section,, 
the unknowns CO,CC^,, 

solving the following system of Nonlinear equations (simul« 
taneously ) ; 


H2 »H20,FgC 12 ,HC1 and NiCl2 may be obtained 


Eqs. (3*12) to (3.15) 
Eqns.(3.16) or (3.17) 
Eqs. (3.18) or (3.19) 
Eq. (3.20) 


The e'.quilibrium 
NiO, Ni, Fe^O^ 
Appendix 3.2. 


and Fe2SiO^ 


ratios at which the solid species 
are stable is given in 
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The Equilibrium constants were obtained from the Free energy 
data, at the temperature of interest. 

Analysis similar to the one delineated in the case of 
Copper Segregation, including the variation of the process 
variables, is employed in the Nickel system too. 

The next chapter outlines the Algorithms employed for 
solving the equations, followed by discussions on the 
results so obtained. 
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CHAPTER 4 

METHODS OF SOLUTION 

As stated in the Introduction, this chapter discusses 2 
Algorithms that generally overcome the deficiencies implicit 
in the solution of a system of simultaneous Nonlinear 
Algebraic equations by the classical methods. 

In the preceding chapter, the problem had been formulated 
as the solution of solving N real equations in N real unknowns; 
the system being given as 

f^iX) ^ fj^(X^,X2,X3, Xj^) =0 

♦ 

♦ 

BROWN’S Method i 

This is a local method^ which handles the 

the 

function (4,1) one at -a time so that £ information obtained from 
working with fj^ can be incorporated when working with f 2 etc. 

A successive substitution scheme is used rather than the 
simultaneous treatment of the f^ characteristic of Newton* s 
method. Brown’s method is derivative free; moreover, second 
order convergence has been proven. 
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Step 1 : 

Let x’^ denote an approximation to the solution X" of 
(4,1). Expand the first function f^ in an approximate 
Taylor series expansion about the point 






(4.2) 


where defined to be the first difference 

y quotient approximation. 


If X*^ is chosen close enough to X'' ,f 3-(X)'^0 ,and solving for 
that variable, say Xj^, whose corresponding approximate partial 
derivative largest in absolute value, we have 




where 
and f^ 


•J 



From (4,3), it may be noted that Xj^ is a linear function of 
the N-1 variables, and may be renamed Lj^( X^^ , X 2 , X^ , X^ , .,,,Xj^^j^), 

Step 2 j 

A new function) 92 of the N-1 variables Xj^,X 2 ,X 2 , 
may be defined, which is related to the second function f 2 f 
as 
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Let g2 be defined as = f2( ^ ^N-1' 

Now expand 92 in an approximate Taylor series expansion about 
the point (x^^x^y ..., linearise (ignoring higher order 

terms) and solve for that variable, say whose correspond- 

ing partial derivative g^y »>,? largest in magnitude s 

N -2 


^N- 


■ 1 “ ^N— 1 ~ ^ ^^2X »h^^2X )* ( x^-X . )-g2/g2)( ) 


where 


o*^ - a (x^ x^+h^ y'^ )-.a'^ 

^2x.jh ~ ^2^^1*^2* ^ ^1+1^ • • * ^N-1^ ^2 

3 - „ 


where is a small positive number. 

Again, = S^-.i(^l »^2 ’•* ^N- 2 ^ 

Step 3 ; 

^3^^1^^2* •••» ~ * * • » ^N— 2^ N— 1’ 

The basic processes of Taylor expansion, linearisation, equating 
to zero and solving for one variable, say Xj^_^2> (whose corres- 
ponding approximate partial derivative g^j^ largest in 

magnitude) is followed. 
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The procedure is continued in this fashion, eliminating 
one variable for each equation treated. Every time a new 
linear expression is obtained, for one of the variables, 

sayXj^^j^, in terms of the remaining N-k-l variables, x^,X 2 »X 2 »*.. 

"this linear expression is used wherever 
had appeared in the previously defined expressions » 

^N-k+2» ^N-1^ ^N* 

Step N ; At this stage, wo have 

^1* ^ 2^ ^ 3^ ^4’’ • * * ” where the Lj’s are obtained 

by back substitution in the N-l rowed triangular linear system 
which now has the form 


L 


x"? - 


i-1 

j=l 


^^N~i+l,x.;h/^N«i+l,x. ;h^ ■ ^ 

J i 


L. 

J 




“ ^N-i+l/^N-i+l,x^;h » ^ - N,N-1, ... 2 

(with gj^ = fj^ and = x^) so that gj^ is just a function of 

the single variable Xj^. Now expanding, linearising and solving 

1’ 


for Xt , we obtain 




Xj^ ? h 


The point x^ thus obtained is used as the next approximation 

11*4* X ^ 

Xj^ , to the first component of the solution vector x • 
Xj^ is renamed as and the system is backsolved to get 
improved approximation to the other components of X , 
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Brown's technique is a fast local technique, that leads 
to convergence if the start values are in a region of local 
convergence^ To check the results obtained from the Brown's 
method, as also to assess the viability of the method for a 
general simultaneous Nonlinear Equation problem, the Hybrid 
method of Powell that places no restriction on the initial 
estimate of the unknowns being 'sufficiently close* to the 
solution of the equations, was used. A general description 
of the method follows. 

POWELL'S Hybrid Method^^®^ ; 

Most Algorithms for the solution of systems of simulta- 
neous Nonlinear equations are centred around either of the 
2 approaches ; 

1. Model expanded as Taylor series and corrections to 
several parameters calculated at each iteration on th-^ 
assumption of local linearity, 

2. Various modifications of steepest descent methods used. 

Both these methods not infrequently fail, Taylor’s series 
method due to divergence of successive iterates, and the 
steepest descent methods because of very slow convergence 
after the first few iterations. Powell's method combines 
both the Newton and steepest descent methods in such a way 
as to give steady progress and a fast rate of ultimate 


convergence 
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To begin the Kth iteration, the estimate of the solution 
x^, a step length and 2 numbers E and M are required. The 

step can be changed on each iteration, and its purpose is ' to 
restrict the length of the displacement ( ) in order 
that the iteration decreases the value of F{x) where F( x) is 
the sum of the square of the Residuals 

N 2 

F(x) = E ifAx)r 
i=l ^ 


E and M determine the conditions for finishing the iterative 
process. It finishes if the value of F(x) is reduced fo 
less than E, or alternatively if the gradient of F{ x) is so 
small that the distance from x to a solution is predicted to 
exceed M. 


First, the kth iteration calculates the elements of the 
Jacobian Matrix 


M _ ^ 




x=x 


(k) 


(4.5) 


and then it evaluates the full Newton-Raphson correction § 
by solving the linear system 


(k) 


f .(x^) + E = 0 


ij “ J 
(k) 


(4.6) 


and also the gradient g 
components 

Si"’ = F(x)]^^^(k) 


,(lc) 


of F(x) at x^ ^ by calculating the 


(4.7) 
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It then tries the test 

F(x^)- > M i i2 (4.8) 

and if it holds, the iterative process ends because of the 
likelihood that the sequence of estimates is converging 

not to a solution of the equations but to a local minimum of 
F(x). 

If the condition (4,8) does not hold, then a displacement 
is calculated to add to the vector x^^^. This displacement 
is just the classical correction if > |( II 2 • 

But if 

< II II 2 (4.8.1) 

then is made equal to 

In this case the displacement has the form 

(4.9) 

where oc^ and are scalars, 

=? o if the step along the steepest descent ‘vector of 
F(x) 

g(lc) = / I (4.10) 

does not go beyond the predicted minimum of F(x) along the 

fk) 

steepest descent vector from x'' , this predicted minimum being 

at the point • 
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xO') - [ ^ / lljf') gC'^lla] g'*'’ (4.11), 

So the condition 

^ l|g'’=^!l^ /l|J<’=^ (4.12) 

2 2 

is tried and if both the inequalities (4,8.1) and (4,12) hold 

then is defined by (4,10). In the case when condition 

(4,8.1) is satisfied, but (4.12) is not satisfied, the point 
f k) -“-f kl 

( x'' + g'' '') is let to be on the straight line joining 

the point (4,11) to the point ( +S^^^) • 

— ('k'} 

Thus ' interpolates between the steepest descent and 
the classical Newton Raphson correction. The functions fj(x) are 
calculated at the incremented value (x^^^ if 


F(x^^^ < F(x^^^) (4,13) 

is satisfied, then ... (4,14) 

and the convergence criterion is tested. However, if (4.13) 
fails, - 


jj(^+l) _ and test (4,15) leads to a reduction 
in the step length 

Revision of the step length A<‘'> constitutes the last 
stage of the iteration which depends on the predicted value 
of the sum of squares of Residuals at (x^^^ 







2 


(4.15) 



which is less than F( x^^'’ ^ ) . If the predicted value is found 
so bad that the actual value of the sum of squares satisfies 
the inequality 

> (l~e) + 6 0^ (4.16) 

where €" is a constant from (0,1), then it is concluded that 
the linear approximations to the functions desired from 

the Jacobian elements (4.5) are not adequate over the distance 
11 §^^^112. Therefore, is reduced by multiplying it 


by a factor ^\, o < A < 1 
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CHAPTER 5 

RESULTS AND DISCUSSION 

In the application of the Algorithms to the Segregation 
problem, the structure of the problem was modified by convert- 
ing the thermodynamic variables (namely the mole numbers) into 
thoir logarithmic values. This was done to prevent one of the 
variables from assuming a negative value (since negative mole 
numbers have no physical significance) and secondly, with a view 
to let the variables to be of similar magnitude in the region 
of interest, and avoid cancellation errors in the evaluation of 
the objective function. 

In the modified form, both Brown’s local method and 

Powell's Hybrid method showed enormous rigidity. Since the 

bane of all methods of solution of Nonlinear equations is the 

initial guess values, these were varied over a wide range of 

—7 

values. Flat start values with as remote a value as 10 upto 
+7 

10 through zero yielded convergence in the case of Powell’s 

method, while Brown's method diverged for very low values of 

—2 —7 +7 

10 through 10 , converging for values from 0 to 10 , Thus, 

between the two, Powell's method showed greater stability. 

This arises from the fact that Powell's Algorithm shares with 

the gradient methods their ability to converge from an initial 

guess value which may be outside the region of convergence of 

other methods. The Algorithm shares with the Taylor series 
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method the ability to close in on the converged values rapidly 
after the- vicinity of the converged values has been reached, thus 

combining the best features of its predecessors while avoiding 
their serious limitations. 

Copper Segregation : 

For the thermodynamic analysis of Cu Segregation that 

require the stoichiometric values desired from experiment, the 

typical values of Yazawa have been made use of. This is 

a lOOg pellet of Cu, consisting of Cu, 205 ^ Fe, 1 /. NaCl, 

l^C (and rest as gangue) . The value of K at ambient pressure, 

and at the temperature of interest is calculated from the free 
(30) 

energy data . The resulting Nonlinear equations were solved by 

the methods outlined in the preceding chapter. The results 

of the computation at a Segregation temperature of llOO^K is 

tabulated in Table 5.1. (for various values of CO 2 /CO ratio), 

maintaining a con'stant value 

The evolved gas consists predominantly of, CO^, H 2 O,/ while the, 
profiles of the other gaseous species vary depending on the 
oxygen potential. In confirm! ty with the very high dissocia- 
tion constant of the monomer - Trimer Chloride Equilibria, the 
partial pressure of the monomer is very low in the Cu Segre- 
gation system. 

The plot (Fig. 5.1) of the partial pressures of the 
gaseous species against the oxygen potential indicates that the 
partial pressure of Cu^Cl^ is high in the region where CU 2 O 
exists, implying thereby that the chlorination reaction (3.5) 
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proceeds forward. It decreases as the oxygen potential 
decreases, while the partial pressure of HCl increases, showing 
the occurrence of eqn. (3.6), yielding copper. The HCl (g) 
thus produced would be recycled and utilised in the chlorina- 
tion of CU 2 O, represented by eqn. (3.5). The gradients of 
the various gaseous species may therefore be directly related 
to the mechanism of the process. It may be noted that among 
the chlorides, Cu^Cl^ and HCl determine the direction of the 
process while FGCI 2 plays no major role. 

The effect of the process variables, namely the amount of 
reductant, water, chloride and the temperature may be studied 
by varying the input compositions of the reactants, and substi- 
tuting the appropriate values of the Equilibrium constants at 
different temperatures. The results of this exercise may be 
plotted as shown in Figs, 5,2 to 5.5, The variations of these 
parameters exhibit a definite trend in the profiles of the 
gaseous species, from which conclusions may be drawn from. the 
thermodynamic standpoint. If an increase in partial pressure 
of Hydrogen chloride over the base value be maintained as the 
yard stick for enhanced Cu segregation, it may be seen that 
decreasing the temperature, decreasing the amount of water, 
increasing the amount of NaCl, and decreasing the amount of 
carbon aids segregation, A higher amount of carbon leads to 
insitu reduction of the oxide, suppressing chloridisation, 
and an excessive amount of water dilutes the effect of the 


chloride 



Logarithmic equilibrium partial pressures of the gaseous species 
at varied log /Prn) values at llOO^K 



rHl>C^'«;tlOC0OCMM3»Has 
C^vOCOCO»HCO\00>‘'<l’COsO 
HCMCOCO '«53''^^\rinvOt^ 

CMcyc^<NCsl CN<N<NCNC'lCN! 

t t ( t t I i I I : I 


oiiOinooiC'JO’^ooja' 

t^’Vt^t^’^O'OCOCOO'CT' 

NO'VrHcpinciCOiDCM 

'^■^'<t(r0r0C0CMW<M>H»-4 

I i i i t t t i 1 t i 


r^O'OO oi>voiOinr-'-J 

OOONCONO»n-^ro(N>H,H 

♦ > ^ ♦ 

I i I I f I I I ( I i 


cotn>~io^coOtoi>-oo^vo 

ro<N<M*H<-ic4CJrO'j30'iO 

'ODOOOqODCOOO/OOCOCOCOO 

OOOOOQOOOOO 

J I I I I T I I I II 

CM>HC4?4ch'Or~C'(N-^vO 

lO-vinrHt^in'VfOcoC'j'-i 

'O'5tf0C0<NOJ(NCM<N<NCM 

% D I 4k ’« •' « e '• • 

OOOOOOQOOOO 

I 1 I T I I T I I I I 


iocoinro>H<MOoiint^aN 
'HinCslCM'^C^’f-tiOO'OOP- 
lO O 03 O Ci. ’'f O O 

I ! T I I I I 1 I I I 


iooo'Viorocor~<-ir-'5j-cM 

cor-OiocNOO'C^cococo 

oot^r^'D'OvOiniOioioio 

OOpOQOOOOOO 

! I T I T I ! I I ! I 


lOoDWin^oob-tHc^'^cM 

oocNjOO(N»f>o^rtocooo 

C0O<N^'Q00OC0^tC0O 

? *HrHfHi fH«-*iC4<NC4C\lOO 

I i I i { I t t I I 


OtnOinOioOinQino 

o lO r- o 04 t> o CM in 

O O O O rH r-f fH CM CM CM 


TJ 

4 -> 

C 

O 

O 



v0 

o 

OJ 

o 

in 

CD 

!>• 

o 

CM 


O 

vO 

in 


in 

vO 

OO 

CO 

CM 


o 

<N 

tn 

00 

H 



r-i 


lO 

• 



'♦ 





u 

» 

CO 

CO 

CO 

(0 




in 

in 

in 

1 

|. 

i 

1 

1 

1 

1 

1 

1 

1 



CO 

r- 

r- 

O 

in 

00 

o 

o 

m 

KJ- 


H 



CM 

o 

o 

o 

vD 

lO 


\]- 

CO 

CO 

CO 

CO 

CO 

CO 


' « 


• 



tt 

♦ 

:h 

♦ 


r-H 

*-1 



rH 


iH 

rH 

«-4 

1 

I. 

1 

i 

1 

1 

f 

1 

1 

1 


CO 

00 

CO 


.H 

o 

CO 

CO 

in 

iD 

\0 

CM 

o 

CO 

c- 

vC 

in 




o 

o 

o 

o 

o 

o 

o 

c> 

CT' 

00 

b 

& 

c 

A 

h 

fe 

«* 

c 

.9 

* 




CO 

00 

CO 

CO 

CO 

CO 

CO 

5 

8 

1 

I 

I 

1 

{ 

8 

5 

5 


CD 

in 

r-1 

in 




CM 

in 

iD 

(M 

rHl 

rM 

t-M 

CM 

CO 

m 

r- 

lO 

in 

O 

iH 

CM 

CO 


in 

vO 

i> 

CO 

00 

h 

9 

• 

« 


« 

« 

% 

• 

ft 

iH 

i—C 

rH 

rH 

rH 

#HI 



rH 

iH 

5 

1 

1 

£ 

J 

1 

1 

1 

! 

1 


00 

H 


oo 


rM 

r^ 

o 


o 

O 

o 

CO 

00 

CO 

r- 


r- 

CM 

CM 

i-M 


r-4 

H 

r-l 

rM 

rH p-M 

ft 

& 

¥ 

'« 

¥ 

♦ 

ft 

♦ 

^ . • 

? 

O 

f 


0 

1 

0 

1 

? 

o 

J 

o 

f 

o o 

8 i 


rH 

CO 


rH 

r^ 

CO 

iH 

00 

r> 


CM 

MO 

o 

in 

o 


O 

00 

€0 

co 

O 

rH 


MO 

00 

rH 

CO 

sO 

00 

rH 

'ft 

¥ 

ft 

¥ 

¥ 

¥ 

ft 

4 

» 

ft 

CM 

3 

CO 

CO 

CO 

CO 





lO 

1 

2 

3 

! 

1 

I 

1 

1 

I 


o 

CO 

i> 

vO 

in 

ID 



CO 


CO 

r- 

o 

r- 

j> 

C*- 

r- 

t-- 

!>• 

r- 

lO 

in 

uo 

in 

in 

lD 

lO 

in 

m 

lO 

■ 3 * 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

• 

o 

o 

o 

CD 

o 

o 

o 

o 

o 

o 

1 

1 

1 



1 

I 

s 

1 

I 


o 

CO 


in 

in 




CO 

CO 

CO 

£> 

CM 


CM 

C-“ 

CM 


CM 

c- 

00 

to 

CO 

o 

CO 

in 

00 

o 

CO 

lO 

ft 

,ft 

ft 

ft 

ft 

* 

ft 

• 

ft 

4 

CO 

CO 

CO 





in 

in 

in 

1 

1 


1 

1 

1 

f 

i 


1 


in o lO o 
r- o (N in 

b « li * 

H| 04 PO CO <»? 


in Q »o © in o 

t- O 


CO 


CM in o- 

h h 

^ ^ ^ 


o 

« 

iD 




f ig.l.! of port lot prmtmrm of th« vorlooi goiooys s^cli« opiimt 

OMygoo f>ot«#it«il of t!00*K fei coppor fogrigcitiofi tyittm. Afiwiifiti 
poftic^ting in t*i« ftoctic»» ort: 

€:a3%» 1^0:0 JV#, H|0: 114% 





‘"S Pco/Pu 


Fi.j-j t .ut of poMia: pressuras against oxygen poteniiol wfwn tr« 
poniepaf.ng amounts of carbon vuneo, other ccmd.f.uf.t, 

%-'mu;nriy tne- sum« cfc tfi Ffg.5,1. 



OM % 
1. 14 % 
2 . U % 








28 


Nick el Segr egation ; 

Calculations analogous to the one employed in the case 
of Cu Segregation may be carried out for the case of Nickel 
too. Here, a 100 g pellet constitutes 2 g Ni, 20 g Fe, 

6 g CaCl 2 j 3g C with the balance as gangue. Incorporating the 
values of Oq , n^^^ and the values of K at the temperature 

1200^K, the equations in the Nickel Segregation system are 
solved simultaneously. The results are tabulated in 
Table 5.2, and a plot of the partial pressures of the various 
gaseous species against the Oxygen potential is made in 
F xg « o « b « 

Contrary to the case of Cu Segregation, the predominant 
species are CO^, HCl and H 2 O, with a fairly significant anjount 
of FeCl 2 too. As in the case of Cu Segregation, the curvp. of 
NiCl 2 exhibits a sharp gradient, but without a consonant/ 
match from the profiles of HCl. The sharp fall in the NiCl 2 
curve indicates that it could be transported to the carbon 
surface to get reducod to metallic Nickel, but the high 
vapor pressure of HCl at the region of higher Oxygen poten- 
tial betrays this conclusion. The chloride vapors may, 
therefore, not be utilised in a cyclic manner for the chlori- 
nation within the pellet, leading to a larger consumption of 
the chloride, a fact that has been corroborated by experimental 
studies by a number of Research workers. 
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The study of the effect of various process parameters is 
made by varying the input composition of the reductant, 
chloride, water, and by varying the temperature. The results 
are plotted in Figs. 5.7 to 5.10. Since the Hydrogen chloride 
generated in the Nickel Segregation system is lost in the bulk 
furnace gases, impeding a cyclic reaction characteristic of 
Cu Segregation, the premise utilised in the Cu Segregation 
cannot be employed for the Nickel Segregation. However, it may 
be observed that the NiCl 2 partial pressure increases with 
increasing CaCl 2 amount, decreasing amount of participating 
carbon and increasing temperature. 

Comparison of the Nickel and Copper Segregation systems on 
a broader schema indicates that Segregation reactions seem 
more favourable in the case of. the latter than the former. 




Table 5.2 

Logarithmic equilibrium partial pressures of the gaseous species at 
varied logCp^p^ /Pnn^ values at 1200°K 
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ANALYSIS THROIGH 


CHAPTER 6 

FREE ENERGY MINIMISATION TECHNIQUES 


In a closed, cheniically reacting system at constant 
pressure and temperature, the potential function which 
governs chemical changes is the Gibbs free energy « The 
condition for chemical equilibrium is that the free energy 
must be a minimum. The free energy minimisation technique 
is based on this principle (32-34,42). 


In the present study, thermodynamic calculation were 
carried out on the system Cu-Fe-C-H-Cl-C using the Free energy 
Minimisation technique. At a given temperature and pressure, 
the free energies of all possible mixtures of a given set of 
species form a surface in the hyper space, the coordinates of 
which are the species, certain areas on this surface cannot be 
realised due to the mass balance constraint not being satis- 
fied, The objective then is to find that location on the 
accessible areas whore the free energy is a minimum. The set of 
coordinates, or mole fractions, that define the m.inimum point 
constitute the desired equilibrium composition. 


Formulation 


The total free energy of a mixture of gaseous and 


condensed species can be expressed as 

n m 

S f ? + I 


-K- 


RT 


i=l 


j=l 


1 


( 6.1a) 
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Thus, thG free encrqv functi 


F(x) = 


on 
n 

S f? + 
i=l ■“ 


j=l 


f! 

3 


(6.1b) 


where , 

X is the set tormeci by the mole numbers of all the species 


Y — f ’^9 '/-Q . './S ,^,,s 

ff = Xf (C| + In , 


\r3 \ 

“ ' y y 


X = 

n 

£ 

X? 


i=l 

1 

II 

O 

F 

rT 

1? + In P 


f ® = X® C . 
J J 3 

P'S F j s 

" RT I j 


For obtaining the equilibrium composition, a set of positive 
mole numbers, satisfying the mass balance and minimising the 
total free-enorgy, are to be generated. At constant temperature 
minimization of free energy is equivalent to minimizing the 
free energy function F. The mass balance constraint can be 
written as 


n 

£ 

i=l 



X 


<3 


m 

+ £ 
j=l 




for e = 1,3, •••, NE 


( 6 . 2 ) 


> i-*-' 



837 . 3 ? 
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where b^, = elemental abundances 
NE = number of elements. 

Let Y be an initial guess for the mole numbers of gaseous and 
condensed species, such as 


Y = (Yf, 


Let Y be a positive set satisfying the mass balance constraint 
The free energy expression for this composition is 

Y? m 

F(Y) - S Y? (Cf + In ) + E y! C® 

where Y = 2 Y? 

i=l ^ 

Let be the difference between X and Y so that 


A? = xf - Yf 
A? =X--Y* 
A = X - V 


By expanding F(X) about Y using Taylor’s series (and neglecting 
higher order terms) 0(X) can be obtained 
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Q^(X) 


1 

2 




1 + E 

dF 

A 

1X=Y 


dxf 

n 

n 

d^F 

E 

E 


i=l k =1 

dxf dx: 



1 

n 

m 


E 

E 

d F 

i=l 

ia 

j=l 

dxfdx? 

m 

£ 

m 

V 

d^F 

j-1 

g=l 

dxjax= 






X=Y 


X=Y J 


A? 

*^1 k 


1 . s .s 

X=Y ^ ^ 


Among the derivatives, 


dp a ^i 

dxf ^ X 


dx| 


= =l 


d^F 


dX^dxg 


X 


1 

X 


1 


if i = k 


if i k 


and 


d^F 


dxfdx? 

^ J 


d^F 


dX®dX® 
3 9 


= 0 


(6.4) 
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Substituting these derivatives in (6.4) 

n Y? 

(2f(X) = F(Y) + £ (cf + In A? 

i=l ^ Y ^ 

+ E ^ 4’.]2 (6.5) 

j=l ^ J ^ i=l ^ Y? Y 

^(x) is to be minimised subject to the side constraint of mass 
balance equations (6.2). Now, another function L(X) may be 
defined combining (6.5) and (6.2) using Lagrarian multipliers 

L(X)=0(X)+ S 

e 1 **“1 

Now to minimise this function, 

»t-(x) ^ aLlx). = 0 
ex? ex! 

i j 

The change of free energy with a change in moles of gaseous 
species becomes 

= [C? + In ^ ] + A - ?] - S p af = 0 (6.7) 

Similarly, for solids, 

= C! - E p a® = 0 

dx= J e 


( 6 . 8 ) 
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Obtaining x| from (6.7) 


X? 


-Yf [Cf + in + y| [L] + a9^] Y? 


(6.9) 


Summing (6.9) over i for all gases. 


E 11 2 a? Y? 
0=1® ^ 


Y? 

E(C? + In -i) Y? 

Y ^ 


( 6 . 10 ) 


n 


Let r , = r, = _ 

ek ke ""le ik 


2: (a?„ a?J Y? 


for e,k = 1,2, ..., NE 

By substituting (6.9) in (6.2), NE number of equations are 
obtained. These equations together with (6.8) and (6,10) 
form a set of (NE+ra+l) linear equations, for the unknowns : 
jig (Lagrarian multipliers) mole numbers of solid species 
X®, X^, X® and the value of X/Y . 

The set of equations to be solved are : 

P^(X/?) + axiX®+a|^X| + ... + a®j_X® + ^INE^^NE ' 

bj, + E a^^f^ 

P2(X/Y) + + ^ 22^2 ^21^^l'*' **• ^2NE^NE “ 

b^ + 


®1NE^1 ••• ®mNE^m ''' ^NEl ^1 *** ^NENE >NE “ 


‘INE'"! ’ •” mNE 


bMG "t 2 
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^NE^E ■" (6«12) 

+ ... + = C® 

s s _ s 

^21^1'^®22^2 -^ ••• + ^2NE^E “ ^2 


®ml^l'^^2^2'^ *•* "^ ^ME^E “ ^m 

^ y9 

where ' 

The solution to the set of equations (6.11), (6,12) and (6.13) 

give a new approximation for the condensed species. The 

values of Pe, ^^e and (x/Y) when substituted into (6,9) yields 

the new values of X?. By repeating this process, the final 

equilibrium value is reached. The method of steepest descent 

was used to facilitate ultimate convergencej this was defined 

to have been reached when the sum of the absolute values of 

the changes in mole fractions between the subsequent itera- 
—5 

tions was < 10 

Analysis of C u2 0"Fe2 0Q“ -C~HCl System 

The input composition of the components chosen was the 
one typically encountered in the segregation process of Oxide 
Copper Ores, From this, the elemental abundances are given by 

^Cu “ ’^Cu^O 

^Fe = ’^Fe203 

b^; = 

^01 ^HCl. 
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where represents the mole fraction of the input species i. 

The gaseous species, the condensed species and the 
elements in the system are to be specified initially. The 
elements in the system are Cu, Fe, C ,H, Cl, 0. The free energy 
data of all the Chemical species conceivable in the system, at 
the temperature and the pressure considered, are to be speci- 
fied. These were obtained from the JANAF tables^^^^* 

A number of gaseous species were considered in the system. 
These include CuCl, Cu^Cl^^ FeCl, FeCl 2 , FeCl^, Fe 2 Ci^, 00, 

CO 2 , H 2 , H 2 O, HCl, O 2 , Cl 2 ^ From the partial pressures 

of the gaseous species so resulting, that have been tabulated 
in Table 6.1, only the following were significant : Cu^Cl^, 

CuCl, FeCl 2 , H 2 » C02» HCl, The condensed species 

that were considered at equilibrium were Cu, C and Fe 203 . 

The following test may be conducted for the presence of 

I 

solid phases. For each of the solids, if the quantity S 

[(F/RT)^ “ ^ ^je^ (6.14) 

is less than zero, then that particular solid is to bo incor- 
porated into the calculation. 

Results and Discussion : 

With the input composition of (in mole fractions) 

CU 2 O = 0.06, Fe 203 = 0.73, C - 0.16 and HCl - 0.05, the final 
equilibrium conposition of the gaseous species and the 



Table 6.1 

Equilibrium mole fractions of the chemical species 
at IlOO^K and input ingredients of CUoO:0,'06,FeoO-:0,73 
C;0.16, HC1;0.05 i. z J f 


So 

SI.' Chemical species Mole fractions 

No. 


1 .- 

CuCl(g) 

2 . 

002013 ( 9 ) 

3 

FeCl(g) 

4. 

FeCl 2 (g) 

5. 

FeCl 3 (g) 

6 . 

Fe^Cl^Cg) 

7. 

CO(g) 

8 . 

OO^ig) 

9. . 


10 . 

HgOCg) 

11 . 

HCl(g) 

12 , 

02 ( g) 

13, 

Cl 2 (g) 

14. 

0 (g) 

15. 

Cl(g) 

16,. 

Cu(s) 

17. 

C(s) 

18. 

Fe^^s ( s ) 


0.1655x10“^ 

0.289xl0’"’^ 

0 . 162 x 10 “^^ 

0^56x10"^ 

0 .-134x10"*^ 

0.634x10“^^ 

0*747x10*“^ 

0.446x10”^ 

0*147x10”^ 

0.887x10“^ 

0.205x10“^ 

0.659x10“^^ 

0.218x10"*^^ 

0.134x10" 

0.12213 

0.765x10**^ 

0.689 
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condensed species were computed by the free energy minimisa- 
tion technique at a temperature range varying from 800*^K to 
1300*^K. The results indicate that a highly reducing atmos- 
phere prevails in the system, and all the Oxide Copper gets 
reduced to Copper. The test for the presence of Cuprous 
well as Cupric oxide failed to yield positive results. The 
test indicated the likelihood of Fe 304 ,Fe 0 and Fe which have to 
be incorporated in future work to check for certain, the 
presence of any of these condensed species. 

A plot of the logarithmic value of the mole fractions of 
the chemical species at equilibrium as a function of tempera- 
ture is shown in Fig. 6.1, Among the gaseous species present 
in significant amounts, CuCl, FeCl 3 and Cu 3 Clg exerted relati- 
vely lower vapor pressures. The vapor pressures of CU 3 CI 3 and 
HCl increased with decrease of temperature. The amount of 
copper extracted predicted from the analysis at equilibrium 
did not show significant perturbation with temperature. although 
from a strictly numerical viewpoint Cu extraction soemed to 
enhance with increasing tomperatuie. Thus while the behaviour 
of HCl with temperature variation is identical to the one 
obtained by Yazawa^^^, the conclusion that an increase in the 
partial pressure of HCl automatically moant enhanced segrega- 
tion is not in confirm! ty with the results obtained here which 
show higher vapour pressures of HCl coupled with lower amounts 
of copper (albeit relative differences in their amounts being 
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very small). Before drawing definite conclusions, the input 
compositions may have to be varied, and all possible condensed 
species included in the analysis to simulate the situation 
in a segregation system. Analysis of the 7 component system 

Cu“Fe-Si-C-H~Cl--0 to take care of the presence of Si02 in the 
Ore would be useful* 
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CHAPTER 7 

LABORATORY SCALE SEGREGATION ROASTING 

Laboratory Scale Experiments on Segregation Roasting 
of synthetic Copper Ore were designed to assess the effect of 
the various process variables on the extraction of copper. 
Synthetic Copper Ore with varying contents of the oxide of 
copper were employed. 

Equipment and Test Procedure ; 

A miniature sketch (not to scale) of the Segregation 

apparatus that was finally used in the laboratory is available 

in Fig. 7.1, the previous attempts at using a 20 mm dia, 75 cm 

and 

long Quartz tube^an Alumina tube having lead to their breakage 
due to 

/exposure to high temperature and moisture. The body of the 
apparatus is of stainless steel, 30 mm dia and 200 mm long. 

Two tubes of 3 mm and 6 mm dia are welded to one of the 
flanges, the former for carrying the moisture, and the latter 
to introduce the thermocouple as also to act as the outlet of 
the outgoing gases, 

A schematic sketch of the set-up for lab scale segregation 
experiments is shown in Fig. 7.2. Argon is purified by- passing 
it through copper gauze maintained at 2G0-300^C in a tube 
furnace, and then allowed to bubble through a water saturator 
at room temperature to provide water vapor in the sample 
atmosphere. 
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In the laboratory scale studies, a modification of the 

# 

experimental practice employed by Opie et was used 

for the treatment of 50g samples of synthetic Ore. In this 

method, which is especially useful due to the small amount of 

synthetic Ore utilised, the copper coated carbon fraction of 
the treated charge is separateci from 

_£the rest of the product by screening on a ASTM 60 mesh sieve. 
The appropriate sizing of the carbon originally, added to the 
reaction charge makes this possible. Initially Segregation on 
Oxide Copper Ore was attempted, but had to be abandoned due 
to the extremely low amount of Cu in the Ore (0 .3/< Cu) that 
lead to difficulties in analysis of the sample after Segrega- 
tion Roasting. Thus a synthetic Ore constituting Cu20(5>i-10;<) , 
Fe 203 ( 10?^-2Q;<) and rest Si 02 was employed for the lab scale 
studies. 

For each test, samples of the synthetic Ore, ground to 
minus 100 mesh were mixed with additions of the requisite 
amount of sodium chloride, and desired amount of coke, sized to 
minus 35 plus 50 meshj the mixture as specified in Table 7.1 
was introduced into the apparatus, and the latter placed in 
the Tube furnace (6.3 cms dia and 44 cms long), maintained at 
the desired temperature { 700-750*^0 ) . Argon bubbled through 
the water saturator provided the moisture, and the Chromel- 
Alumel thermocouple introduced through the 6 mm dia tube aidod 
in recording the temperature. The sample was Segregation 
roasted for 1 to 5 hours, at the end of which the apparatus was, 
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Representative results of the laboratory scale Segregation 

experiments 
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removed from the furnace, cooled to room temperature under 
the atmosphere gas, before being removed for chemical 
analysis* 

Results ; 

Screening of the treated charge on a 60 mesh sieve 
resulted in separation of segregate, composed of coke 
particles coated with copper metal, from a fine oxide 
residue. The segregate is made of relatively spherical 
particles of coke coated with spongy deposits of metallic 
copper. The plus 60 mesh portion was subjected to X-ray 
analysis to check for the presence of copper* 

Under the conditions specified, X-ray analysis indicated 
the presence of pure metallic copper. Examination of the 
residual - 60 mesh portion too to X-ray analysis showed that 
small amount of spongy copper broke off and passed through 
the sieve. Such fine particles would be generally recovered 
as part of the segregate when flotation-separation is employed. 

The best results in the lab scale segregation experiments 
were obtained with a mixture consisting of lyi NaCl, coke 
and bubbling of Argon through water saturator at the rate 
of 60 drops per minute. 

Further work, as planned during the start of the present 
study, could not be pursued owing to the snag involved in a 
thorough quantitative analysis of the segregate. 
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CHAPTER 8 

CONCLUSIONS AND SUGGESTIONS FOR FUTURE WORK 

Theoretical and experimental studies in the present work 

showed that : 

1) Qualitatively speaking, Segregation is more favourable 
in the case of copper than Nickel, 

2 ) Brown's and Powell's methods when employed even with poor 
initial estimates lead to the solution of the system of 
the simultaneous nonlinear algebraic equations, 
characterising the segregation system* 

3) From a thermodynamic standpoint, copper segregation is 
enhanced by decreasing the temperature, decreasing the 
amount of water and increasing the amount of NaCl, 

The analysis reveals that the partial pressure of NiCl 2 
increases with increasing temperature, increasing CaCl 2 
amount, and decreasing the amount of carbon. 

4) Synthetic Ores may be employed in laboratory scale experi- 
ments for investigating the effect of the process varia- 
bles in the segregation process. Synthetic Cu Ore used in 
the present study responded to segregation (with the 
modified sizing analysis employed.) 
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5) Application of a generalised computational scheme 
through the minimisation of Gibbs Energy to the Cu 
segregation system yielded the equilibrium composition 
of the chemical species. The teehnique is particularly 
useful due to its inherent ability to calculate the 
equilibrium amounis of the condensed species. This is 
utilised to predict the influence of various process 
parameters on the segregation of Copper . 

Suggestion for Future work ; 

1. The algorithms of Brown and Powell may be applied for a 

thermodynamic study of the systems involving complex chemical 

equilibria , A check on the rigidity of the algorithms may be 
undertaken by 

/employing a varied range of initial estimates of the unknowns. 

2. The effect of a number of process parameters including 
temperature, nature of the reductant, salt, presence of 
moisture* silica, time. Ore impurities, etc, may be studied by 
employing a synthetic Ore to obtain an insight into the mecha- 
nism of the process. 

3. The equilibrium mole fractions of the chemical species 
calculated by the minimisation of Free energy may be compared 
with the actual experimental results to obtain the efficiency 
of extraction of the segregated metal. 
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4. A number of process parameters may be varied^ and its 
effect at the final equilibrium situation may be ascertained 
theoretically by the Free energy minimisation technique* As 
stated earlier, this could be compared with the experimental 
runs so that a juxtaposition of the two would lead to a 
better understanding of the system under differing conditions, 

5. The condensed species that are likely to be present at 
equilibrium as per the formula (6,14) in the Cu segregation 
system for the input composition of the reactants considered 
in our study were Fe^O^, FeO and Fe« These have to be incor- 
porated and the analysis carried out. Inclusion of Si02 in 
the input ingredients would be ideal although the equilibrium 
solution to the problem would then be much more complex, due 
to the probability of the presence of a number of solid 
specie*, 

6. A procedure similar to the one employed in the case of 
Copper may be tried for other potential metals amenable to 
Segregation Roasting like Antimony, Bismuth, Cobalt, Gold, 
Lead, Nickel, Palladium, Silver and Tin. 



The equilibrium constants of the various reactions are listed below 
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Equilibrium Log P 


CO2 



at 


different temperatures 


Chemical species 
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1100°K 

1200°K 
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" f ^ 

1, CU 2 O/CU 
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2, Fe 2 Si 0 ^/F 6304 
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2. 15 

3 . Fe 90 ^/Fe 203 

4,53 

4,378 
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4. NiO/Ni 

2.58 
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2.260 

2.13 
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PROGRAM^Tp SOLVE A SET OF HIGHLY NONLINEAR EQUATIONS: 
USING BROWNES TECHNIQUES 

X( 8 ),P( 8 ),QC 8 ) 

COMMON RR 
COMMON C1,C2,C3 
READ(20j.*)Cl,C2,C3 
ITEMPS1200 

WRITE(21,ll) 

WR|TE(2142),cl,c2,c3atemp 

WRITE(21,13) 

WRITE(21,17) 

WRITE(22,14j; 

WRITE( 22 , 15 Lc 1 ,C 2 ,C 3 .ITEMP 

WRITE(22,165 

FORMATC2C/),T6,?REMARKSS . Varying C ^ 

format(2(/3,t6, 'WHEN Cl»»,F8.5,2X, 'C2»* »F8.5,2X# 
1'C3=%F8.5,2X,^H1 SEGN TEMPt = SI5) 
FORMAT(T6r65(*-')3 

rORMAT(2(/3,T16, 'REMARKS; - VARYING C2') 
FORMATC2C/).Ti6,'WHEN C1«1,F8.5,2X, 'C25' ,F8,5 , 

12X, 'C3s',F8.542X,'Ni SEGRN TEMPS=',I5) 
F0RMAT(T16,65(!-'3) . . , . . ^ 

FORMAT ( // f T4 64 'COS , T31 , 'C02 ' , T46 , ' H2 ' , T 6 1 , ' H20 ' , T76 , 

1 'HCli',T9{,'NlC12',T106,VeC12*,/) 

DO 2 RRsO. 0,4. 0,0.25 . , 

OPENCUNlTaSS, DEVICE® 'DSK', FILES 'FOR35.DAT') 

REA0(2O,») (XCI),lal, 8 > 

TYPE*,X 

DO 10 lsl ,8 

XClisALOGlOCXCD) 

CONTINUE 

TYPE^.X 

CALL NONLINC8, 8,100, i,X,l.OE-08) 

TYPE*,X 

CALL EVALCX,P,Q) 

CLOSECUNITsSs) 

OPEN (UNIT»20,DEVICEs'DSK*,FILEs'F0R20.DAT'3 

WRITE(20,»)X 

CL0SE(UNtT=?20) 

CONTINUE 

STOP 

END 

N NUMBER OF E0UATIONS(»NO, OF UNKNOWNS). 

SUMSIG « NO, OF SIGNIFICANT DIGITS DESIRED. 

HAXIT ^ MAXIMUM. NO.- OF ITERATIONS TO BE ALLOWED. 
IPRINT « OUTPUT OPTION, OUTPUT IF s If 
X .. VECTOR OF INITIAL GUESSES, 

EPS - CONVERGENCE CRITERION. 

SUBROUTINE NONLIN(N,NUMSlG,MAXlT.IPRINT,X,EPS) 
DIMENSION X.C30),PART(30),TEMP(30),COEC30,3U ,ISUfl(30) 
DIMENSION LOOKUP(30,30),GNDREC5) 
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04400 
t4500 
04600 
04700 
04800 
04900 
05 0 00 
05100 
05200 
05300 
05400 
05500 
05600 
05 7 0 0 
05800 
05900 
16000 
16100 
6200 
16300 
16400 
6500 
6600 
6700 
06600 
6900 
7000 
17100 
17200 
17300 
(7400 
17500 
17600 
17700 
17800 
7900 
8000 
8100 
«200 
8300 
8400 
8500 
8600 
8700 
8800 
8900 
9000 
9100 
9200 
9300 
9400 
9500 
9600 
9700 
9600 
'9900 

moo 

OlOO 

10200 

10300 

^400 


IW 


DELTA WILL BE A FUNCTION OF THE MACHINE AND THE PRECISION USED 
DELTAS l,E-07 

RELC0N=10,0E+0»»C-NUMSir*) 

iJTESTsl 

imPRINT.EO.4 3 WRITE(35,48) 

FORMAT (1 HI) 

DO 700 M=1,MAXIT 

I0UIT=0 

FMAXsO. 

M1=M-1. 

IFdPRINT.NE.l) GO TO 9 

FORMAT (//T41 ,157 (T41,6E18.83) 

DO 10 J!sl,N 
LOOKUP C I, J 3**3 

THE ARRAY LOOKUP PERMITS A PARTIAL PIVOTING EFFECT WITHOUT 
HAVING TO PHYSICALLY INTERCHANGE ROWS OR COLUMNS, 

DO 500 K=l,N 
1F(K-13 134,134,131 
KMIN*K-1 

CALL BACKCKMIN,N,X,ISUB, COE, LOOKUP) 

SET UP PARTIAL DERIVATIVES OF KTH FUNCTION. 

CALL AUXFUNCX,F,K) 

FM AX* AMAXl ( FMAX , ABSXF 3 3 

IP(ABSCF).GE,EPS3 go to 1345 

IQUIT*IQUIT+i 

IPCIQUIT.NE.N) GO TO 1345 

GO TO 725 

FACT0R*0.001E+O0 

ITALLY*0 

DO 200 I*K,N 

ITEMP*LOOKUP(K,I) 

HOLDaXCITEMP) 

PREC*S,0E-6- 

PFEC IS A FUNCTION OF THE. MACHINE SIGNIFICANCE, SIG. AND SHOULD 
BE COMPUTED AS PRECaS.^lO. ♦4!(-siG+2) . IN THIS INSTANCE WE WERE 
DEALING WITH AN 8 DIGIT MACHINE. 

BTAaFACTOR*ABSCHOLD3 
HsAMINl (FMAX, ETA) 

IFCH.LT.PREC) H=PREC 
X(ITEMP)*HOLD+H 
IFCK-13 161,161,151 

CALL BACK(KMIN,N,X,ISUB, COE, LOOKUP) 

CALL AUXFUN(X,FPLUS.K) 

PARTCITEMP)a(FPLUS-F)/H 

X{ITEMP3*HOLD 

irc ABSCPARTCITEMP)).LT.DELTA ) go T0.190 
IFC A8SCF/PARTCITEMP)).LE,t.0E+15 3 GO TO 200 
ITALLY=ITALLY+1 
CONTINUE 

IF(ITALLY.LE,N-K3 go to 202 
FACTOR*FACTOR*10 . 0E4-00 
IFCFACTOR.GT.II.) GO TO 775 
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10500 

tioeoo 

10700 

10800 

10900 

11000 

11100 

11200 

11300 

11400 

11500 
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12000 

12100 

12200 

12300 

12400 

12500 

12600 

12700 

12800 

12900 

13000 

13100 

13200 

13300 

13400 

4^3500 

13600 


202 


203 


209 

210 


C 

C 


13700 

13800 

220 

13900 

500 

14000 

14100 

C 

14200 

14300 

14400 

14500 

14600 

610 

14700 

14800 

C 

14900 

15000 

625 

15100 

15200 

630 

15300 

15400 

15500 

649 

15600 

650 

15700 

660 

15800 

700 

15900 

16000 

1753 

16100 

16200 

16300 

1763 

46400 

H^SOO 



GO TO 135 

IFCK.LT.M) GO TO 203 

IF( ABSCPARTCITEMP)),tT, DELTA ) GO TO 775 
cnECK,N+n=o.OE+oo 
KMAXbITEMP 
GO TO 500 

FIND PARTIAL DERIVATIVE OF LARGEST ABSOLUTE VALUE 

KMAX=LOOKUPCK,K) 

DERMAX=ABSCPART(KMAX) ) 

KPLUS=K+1 
DO 210 isKPLUS^N 
JSUBsLQOKUP(K,I) 

TESTS ABS ( PART (J5UB)) . 

IF(TEST,LT.DERMAX) GO TO 209 
DERMAX=TEST 
LOOKUP ( KPLUS , I ) sKMAX 
KMAXsJSUB 
GO TO 210 

LOOKUP C KPLUS , I ) = JSUB 
CONTINUE 

IF( ABS(PARTCKMAX)).EQ.0.0) GO TO 775 

SET UP COEFFICIENTS FOR KTH ROW OF TRIANGULAR LINEAR SYSTEM USE 
TO BACK - SOLVE FOR THE FIRST K VALUES OF Xfl) 

ISUBCK)=KMAX 
COECK.N+DsO.OE+OO 
DO 220 JsKPLUS.N 
JSUBsLOOKUp ( KPLUS . J ) 

C0E(KWSU8)s-PART(JSUB)/PARTCKMAX) 

COECK,N-H)sCOECK,N+l)+PART(JSUB)»X(JSUB) 

CONTINUE 

C0ECK,N+1>s(C 0ECK,N41>-F3/PART(KMAX)+XCKMAX) i 

BACK SUBSTITUTE TO OBTAIN NEXT APPROXIMATION TO X 

X(KMAX)sCDE(N,N*1) I 

IFCN.EQ.l) GO TO 610 

CALL BACK (N-lfNfX.ISUB, COE, LOOKUP) 

IF(M-l) 650,650,625 

TEST FOR CONVERGENCE 

DO 630 I=?1,N 

IF( AB5CTEMP(I)-X(I)>.GT.ABSCX(I))»RELC0N ) GO TO 649 

CONTINUE 

JTEST=JTEST+1 

IF(JTEST-3) 650,725,725 

JTESTsl i 

DO 660 Isl,N I 

TEMP(I)sXCl) 

CONTINUE 

FORMAT (/f4nfN0 CONVERGENCE, MAXIMUM NUMBER OF ITERATIONS USED'! 

IFCIPRINT.NE.t) GO TO 800 ! 

WRITE(35,1763) . . [ 

F0RMAT(/T41,1fUNCTI0N VALUES AT THE LAST APPROXIMATION FOLLOW'/; 

GO TO 7777 1 



4 


46600 
J16700 
16800 
16900 
17000 
17100 
17200 
17300 
17400 
17500 
17600 
17700 
17800 
17900 
18000 
18100 
I 18200 

18300 
18400 
18500 
18600 
18700 
18800 
18900 
19000 
19100 
19200 
19300 
19400 
19500 
: 49600 
^ 19700 
19800 
^ 19900 
20000 
20100 
20200 
20300 
20400 
20500 
20600 
20700 
20800 
20900 
20950 
21000 
21100 
21125 
21137 
21150 
21200 
21250 
21262 
21268 
21271 
21273 
21275 
21300 
21500 
21600 
'fl700 


725 

7777 


7788 


8777 

751 


7515 


7525 

800 


IFCIPRINT.NE.I) go fO 800 

DO 750 K=1,N 

CALli RllXFlJN(X,PARf(K).K) 

CONfINUE - 

IFCIFLAG.NE.IJ GO TO 8777 

FORMATC (T41 ,6E20,8)) 

60 TO 800 
WRITEC35,751) 

F0RMAT(//T41, 'CONVERGENCE HAS BEEN ACHIEVED. THE FUNCTION VADUES 

WRITEC35.7515) (PARTCK) ,Kal ,N3 . .. 

rORpTCTll, 'AT THE FINAL APPROXIMATION FOLLOW; '//CT41 .6E20, 8) ) 

GO TO 800 

WRITE(35,752) 

IS SINGULAR. TRY A DIFFERENT* 3 

WIaI TE C 3 S f T 5 / 

F0RMATCT41, 'INITIAL APPROXIMATION.') 

MAXITsMl+1 

RETURN 

END 


SUBROUTINE BACK(KMIN,N,X,ISUB, COE, LOOKUP) 

THIS SUBROUTINE BACK-SOLVES THE FIRST KMIN ROWS OF A TRIANGULAR^ 
LINEAR^SYSTEM FOR- IMPROVED X VALUES IN TERMS PREVIOUS ONES. 
DIMENSION XC30) ,COEC30,31),ISUBC30),LOOKUP(30,30) 

00 200 KK*1,KMIN 

KM*KMIN-KK+2 

KMAXxISUBCKM-I) 

X(KMAX)=0.0E+00 
DO 100 J=KM,N 
JSUBsLOOKUP f KM , J ) . 

X(KMAX)»X(KMAx5+C0ECKM-l,JSUB)»XCasUB) 

CONTINUE . 

X(KMAX)xXCKMAX)+COEfKM-l,N+l) 

CONTINUE 

RETURN 

END 


SUBROUTINE AUXFUNCX, t,K) 
DIMENSION XC8) 

COMMON RR 
COMMON C1,C2>C3 


data all,AL2/2.l3,2,l5/ 


DATA ALl,AL2/2,26.2,287/ 
DATA AL1,AL2/2,41,2.45/ 


AKla.116; AK2«2.38E3f AK3s4,53? 

AK4«,267; AK5s.§94, 

AK2S1.48E4; AK3a.l74 
AK4=0,0877j AK5S1.09 

AKi=4,7E-2; AK254.96E3; AK3al.07? 

AK4a0.163? AK5=0.)83 


AK3a.l74? 


GK2aALOG10(AK2); 

GKiaALOGlOCAKl); 

GK5aALOG10<AK5)? 


GK4aALOG10CAK4) 

GK3aAL0G10CAK3) 


s 


21800 

721900 

22000 

22100 

22200 

22300 

22400 

22500 

22600 

22700 

22800 

22900 

23000 

23100 

23200 

23300 

23400 

23500 

23600 

23700 

23800 

23850 

23900 

23950 

23990 

23995 

24100 

24112 

24123 

24125 

J4150 

24162 

25000 

25100 

25200 

25300 

25400 

25500 

25600 

25700 

25800 

25900 

26000 

26100 

26200 

26300 

26400 

34300 

34400 

34500 

34600 

34650 

34675 

34700 

34800 

34900 

35100 

35150 

35175 

35200 


2 

3^ 

4 

5 

6 

7 

10 

31 

8 


21 

77 


22 


11 


3 

4 

5 


GO TO Cl ,2, 3, 4, 5, 6,7,8), K 

Y“CCI-(10.*X(1)+10.*X(2))) 

RETURN 

YsCC2-C10.*XC4)+l0.#X(3)+0,5»10.*XC5))) 

RETURN 

Y=CC3-C0.5*10.''X(5) + IO.^XC7) + 10.^XC6))) 

RETURN 

Y=(C10.*XC8))'>C10.fXCl)’l-10.»X(2)+lO.SX(3) + 10.‘ 
l + 10,«XC5)+l0,*X(6) + l0,/^X(f ))) 

RETURN 

Xs^?.S:^^^2)-xci))> 

RETURN 

Ya!CGK5-(XC2)+XC3)-X(l>-X(4))) 

RETURN 

IFCRR.GE.ALl) GO TO 10 
Y»(GK2-(2.*X(5)-X(6)-X(3))) 

6^ To 3 1 

YiCGKl-txC6)+X(4)-2.*X(5))) 

RETURN 

IF (RR,l.T,AL2)G0 TO 73 , 

Y=(GK3-(3.*X(7)+4.fXC4)-6.*XC5)-X(3))) 

GO TO 77 

YaCGK4-(XC7)+X(4)«2.»X(5))) 

RETURN 

END 

SUBROUTINE EVAL(X,P.O) 

THIS SUBROUTINE EVALUATES EQBM S STOICH CONST 
DIMENSION X(8),P(8),QC8) 

COMMON RR 
DO 22 I»l,8 
XCDslG.-XCl) 

CONTINUE 

DO 11 1=1,7 
P(I)=X(I)/XC8) ^ 

■(I)=ALOG10CPCI)) 

ONTINUE 


WRITEC21,l)RR,CQ(I),Isi,7> 
WRITEC22,2)RR,CX(l5,Ifl,45 
WRITEC22,3)(X(I),I=i,8) 
WRITE(22 ,4) CPC I) ,1=1,4) 
WRlTEC22,5)CPCl),ia5,7) 

F0RMATC/,T6,F4.2,7F15.6) 

F0RMATC//,T16,F4.2,T20,4E20,8) 

FORMATCT20,4E20.8) 

FORMAT C /,T20,4E20 . 8) 
F0RMATCT20,3E20,8) 

END 


XC4) 
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00001 

c 

PROGRAM FOR THE SOLUTION OF A SYSTEM OF NONLINEAR 

Jooil 

00021 

c.. 

SIMULTANEOUS ALGEBRAIC EQUATIONS BY THE HYBRID 

c ♦ ® 

METHOD DUE TO POWELL 

ooo:ii 


C1,C2,C3 SSTOICHIOMETRIC VALUES 

RASLOG P'C02/P'CO 

00041 

c ® ® 

00051 

On® 

00100 


dimension X(9),F(9), AJINV(9, 9), W(250) 

OOlOl 


COMMON RA 

00141 


DATA Cl, C2, 03/0.025,0,0633, 0,0154/ 

N=9f MAXFUN=900; a5cs1,E-08 

00181 


00221 


STEPsSQRTCFTOL) 

00261 

00381 


0MAX=5,0 

00421 


OPEN(UNIT=30,DEVICEa'DSK* , FILES *E,CDR') 

00431 

c « » 

FEED THE FLAT START VALUES 

00441 

c » « 

AND THE RATIO OF LOG P^C02/P'C0 

00461 


READ(30,4) CXCI),I=1,8) 

00501 . 


READU0,») RA 

00511 

Cp* 

X(9) IS TOTAL NO OF MOLES 

xb)=o. 

00541 


00581 


DO 5 1=1,8 

X(9)=X{95+XCI) 

00621 


00661 

00671 

5 

CONTINUE 

00701 


DO 10 1=1,9 

X(1)=ALOG10(X(I)) 

00741 


00821 

00900 

iO 

CONTINUE 


CALL NONLIN C 2 , X , F , AJ IN V , STEP , DMAX , ACC , MAXFUN , 1 , W ) 
CALL EVAL(X,P,Q) 

CLOSECUNIT=30) 

00910 


1 00920 


41000 


CALL EXIT 

1 011 00 

1 01200 


END 

! 01300 


SUBROUTINE GALFUN(N,X,F) 

01400 


DIMENSION X(9),F(9) 

01401 


DIMENSION XC9),RCC9) 

01403 


COMMON PA 

01404 

c® • 

AL'S REPRESENT STABILITY BOUNDARIES 

0 1 4 0 5 

c 

DATA Cl, C2, 03/0,0307,0.071 1,0, 0154/ 

DATA Cl, C2, €3/0,025,0.0633, 0.0154/ 

DATA ALl,AL2,AL3/4.67,2,45,4,38/ 

01407 


01409 

01411 

c 

01412 

C »'» 

AK»S REPRESENT EQBM CONSTANTS 

01413 


AK2=92,5; AK4=,174J AK6=1,09; AK7=3,43E+10 

AK1=467,; AK3s8.72,* AK5=,0877 

GK2=AL0G10(AK21? GK4=AL0G10CAK4) 

01415 


01417 


01419 


GK6=AL0G10(AK6)? GK7=AL0G10 (AK7 ) 

01421 


GKl =AL0G1 0 ( AKl 3 ? GK 3=ALQG 1 0 C AK3 ) 

01423 

01425 

01429 


GKBSALOGIOCAKS); 

01430 

c • • 

F'S REPRESENT THE FUNCTIONS WHOSE ZERO'S ARE REQD 
FC13=CI-(10,''(X(1)) + 10.'*XC2)) 

01431 


01433 


F(2)=C2-(10,*X(3) + 10.!^X(4)+0.5*10.'‘XC5)) 

01435 


01437 


P(3)=C3-(10.'*X(5)+2,0*10.*X(8)+3.0»10.?XC7) + 10,^XC63 

01439 

01441 

01443 



FC4}a(10.*X(9))-(10.'fX(l)tl0.fX(2) + 10,*XC3 3 + l0.»X(4) 
l-H0.fXC5) + 10,*X(6) + J0.^XC7) + 10,'‘XC8)) 

01445 

V 



^1447 



©1449 


F(53=KA-(XC2)-XCl)) 

' 01451 

1 01453 


F(6)=GK2-(2,*«'X(5)-(2./3,)*XC7)-X(3)-Cl ,/3,)»XC9)) 

i 01455 

C 

FC6)=GK1-CC§./3.)*X(^)+X(43+C1./3.)*X{9)-§.*XC5)) 

i 01457 


F(7)=GK4-C3,*XC8)+4.4X(4)-6.*XC5)-X(333 

Ft7)=GK5-(xj8)+XC4)-2,4X(5)3 

! 01459 

C 

! 01461 


j 01463 

C 

F(7 3=GK3-(2,»X(8)-i-3.>I'^(4)-4,=I'X(5)-XC3)) 

01465 

01467 


F(8)=GK6-(XC2)+X(3)-X(1)-X(4)) 

, 01469 
! 01471 


F(9)sGK7-(xC7)+2.*XC9)-3.*X(6)) 

01473 


RETURN 

01475 


END 

01477 

01479 


SUBROUTINE EVAL(X,P,Q) 

01480 

C • • 

THIS SUBROUTINE CALCULATES THE EQBM AND STOICHIO VALUES 

01481 

01483 


DIMENSION XC9),P(9),Q(9) 

01485 


COMMON RA 

01487 

c 

DATA ALUAL2,AL3/4.67,2,45,4.378/ 

ITEMP=1100 

01489 


01491 


DO 22 1=1,9 

XCDslO.-'Xd) 

01493 


01495 

22 

CONTINUE 

01497 

01499 


DO 11 1=1,8 

01501 


P(I3=X(I)/X(9) 

01503 


Q(I)=ALOG10CPCI)) 

41505 

^1507 

11 

CONTINUE 

01509 


WRITE(20,7) 

01511 


WRITE(20,83 

01513 


WRITEC20,4),XCn,P(l),Q(l) 

01515 


WRITEC20,51,X(2),PC2),0C2) 

WRITE(20,6);x(3),Pn),0(33 

WRITEC20,15)X(4),PC4),0(4) 

01517 


01519 


01521 


WRITEC20,2)XC5),P(5),6C5) 

WRITEC20,35) XC6) ,PC65 ,QC63 

01523 


01525 


NRITE(20,3)X(7).P(7),0(7) 

WRITEC20,25),X(8),PC8),0C8) 

01527 


01529 

01531 

C 

IF (RR.GE.ALl)GO TO 71 

01533 


F2=X(5)*2/(XCb^C2./3.)»X(3)*X(9)i*(l,/3,3) 

01535 


TYPE*,F2 

01537 


WRITE(20,55),F2 

01541 

01543 

C 

Fl = (XC7)''(2,/3,)*XC4)»X(9)?€i./3.)3/XCS)«2 

01545 

C 

TYPE »,F1 

WRITE(20,53),P1 

01547 

c 

01551 

01553 

c 

IF CRB.LT.AL2)GO TO 73 

01555 

c 

IF (RR.GE,AL3)G0 TO 74 

01557 



01559 

c 

F4=CX(8)'*3#XC4)*4)/£X(5)*6»X(3)) 

01561 

c 

TYPE»,F4 

01563 

c 

WRITEC20,65) ,P4 

01565 

c 

GO TO .77 

, 11567 
%1571 


F5=XC8)4X(4)/CX(5)*2) 
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i)1573 

^1575 

01579 

01583 

01585 

01587 

01591 

01593 

01595 

01597 

01599 

01601 

01603 

01605 

01607 

01609 

01611 

01613 

01615 

01617 

01619 

01621 

01623 

01625 

01627 

01629 

01631 

01633 

01635 

01637 

41639 

01641 

01643 

01645 

01647 

01649 

01651 

01653 

01655 

01657 

01659 

01661 

01663 

01665 

01667 

01669 

01671 

01673 

01675 

01677 

01679 

01681 

01683 

01685 

01687 

01689 

01691 

01693 

01695 

^1697 

^900 


C 

c 

c 


92 

88 

87 

86 

85 

75 

55 

65 

l| 

45 

46 
13 
7 

8^ 

25 

15 

2 

3 

4 

5 

6 


TyPE=»',F5 

'WRlTE(20r6l),F5 

F3-CX(8)*2tX(43*3)/(X(5)*4*X(3)) 

TyPE»,F3 

WRITE(20,62)»F3 

F6s(XC2)*X(3))/(X(UfX(4)) 

TYPE*,F6 

WRITE(20,75)rF6 

F9s(XC7)tX(9 3‘‘2 3/X(6)^3 

TyPE*,F9 

WRITE(20,903,F9 

F7a:X(2)/X(l) 

WRn’E(20,85),F7 

TypE*,F7 

CC1=XC1>+XC2)^ ^ 

WRITE(20,86) ,CC1 

CC2=X(3)+X(4)+0,5»X(5) 

WRITE(20,87) ,CC2 

CC3=X(5)+2,»XC8)+3.»X(7)+XC6) 

WRITE(20,88),CC3 

WRITE(20,45),XC93 
WRITE(20,46),RA 
WRITE (20,1 3 ),1TEMP 

RETURN 

F0RMAT(/,25X, ^MONO-TRIMER EOM CONSTANT!?' ,1PE20, 8) 
F0RMAT(/,25X,'SUM PROM Cl BALANCES? FI 3 ,8) 
F0RMAT(/,25X,*SUM from H2 BALANCES?' , FI 3, 8) 
F0RMAT(///,25X,'SUM FROM C BALANCE??* ,F14. 8) 
F0RMAT(/,25X,'PC02/P CO ratios', F13. 5) . 

FORMAT (/,25X,'C0-H2a/C02-H2 EQ8M CONSTAT*' , FI 3, 8) 
F0RMAT(///,25X, *Co3Cl3 REDN EOM CONSTs' ,F15.85 
FQRMATC///,25X, 'CU OXIDE CHLN EO CONS=' ,F14.8) 



JC A V A # A 4. VOWI9 £.v ^ ^ f ' 

FORMATC///,25X, 'TOTAL MOLESCCALO*' ,E19.8) 
FORMAT(/,25X,'LOG P C02/P CO RAT= ,F8,|3.^ 
FORMAT(/,25X,'SEGREGN TEMPERATURES' .17, 
FORMAT(10(/),T42, 'No OF MOLES ' ,T66 , *P PRESSURES' .T90 , ' 
FORMAT(T42,'————-—',T66,*-———--—'.T90,>— 
F0RMATC//,25X,'CuCl',T4l,El5,8,T65,El5,8,T90,Fll,8) 
FORMAT(//,25X, 'FeC12*,T41,E15.8,T65,ElS.8,T90.Fli.8) 
FORMATC//,25X,'H2O',T4l,El5,8,T65,El5.8,T90,Fil,8} 
FORMAT(//,25X, 'HCL'|T41,E15,8,T65,E15,|,T90,F11.8) 
FORMAT(//,25X,'Cu3ci3',T41,El5.8,T65,E15.8,T90,Fll.8) 
FORMATC//,25X,'CO^.T41,E15,8,T65,E15,8,T90,F11,85 
F0RMATC//,25X.'C02*,T4i ,E15.8,T65.E15.8,T90,Fli.|) 
F0RMATC//,25X. 'H2',T4l,E15,8,T65,E15,8,T90,F11.8) 
CLOSECUNITsig) 

END 


LOG P P|i 
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2000 

2100 

02200 

02300 

02400 

02500 

02600 

02700 

02800 

02900 

03000 

03100 

03200 

03300 

03400 

03500 

03600 

03700 

03800 

03900 

04000 

04100 

04200 

04300 

04400 

04500 

04600 

04700 

04800 

05000 

^100 

05200 

05300 

05400 

05500 

05600 

05700 

05800 

05900 

06000 

06100 

06200 

06300 

06400 

06500 

06600 

06700 

06800 

06900 

07000 

07100 

07200 

07300 

07400 

07500 

07600 

07700 

07800 

07900 


C 

C 

c 

c 


c 

c 


c 

c 


85 

86 


8100 


3 

6 

7 


8 

9 

5 

t 

’ 4 
10 
20 
12 

14 

16 

17 


SUBROUTINE NONLINCN, XtF, AU1NV,DSTEP,DMAX, ACC ,MAXFUN, IPRINT ,W) 
DIMENSION XC9),F(9)„AJINV(N,N),W(9) 

SET VARIOUS PARAMETERS 
MAXC =0 

'MAXC' COUNTS THE NUMBER OF CALLS OF CALFUN 
NT = N+4 
NTEST — NT 

'NT* AND 'NTEST' CAUSE AN ERROR RETURN IF FCX) DOES NOT DECREASE 
DTEST = FL0AT{N+N)-0,5 

'DTEST' IS USED TO MAINTAIN LINEAR INDEPENDENCE 

NX = N»N 

NF = NX+N 

NW a NF+N 

MW a WW+N 

NDCs MW+N 

ND a NDC+N 

THESE parameters SEPARATE THE WORKING SPACE ARRAY W 
FMIN a 0, ^ 

USUALLY ’fmIN* IS THE LEAST CALCULATED VALUE OF FCX3# 

AND THE BEST X IS IN W(NX+1) TO W(NX+N) 

DD a 0, 

USUALLY 'DD' IS THE SQUARE OF THE CURRENT STEP LENGTH 
DSS = DSTEP*DSTEP 
DM a DMAX^PMAX 
DMM a 4,*DM 

^fs'”C0NTRQLS A 'GOTO' STATEMENT FOLLOWING A CALL OF CALFUN 
TINC a 1 

'TINC' IS USED IN THE CRITERION TO INCREASE THE STEP LENGTH 
IF (IPRINT)1,1,85 
PRINT 86 
FORMAT Cl HI) 

CALL THE subroutine CALFUN 

MAXC a MAXC+1 

CALL CALFUNCN,X,F) 

test for CONVERGENCE 

FSQ a 0 

DO 2 lal.N 

FSQ a FS0+FCI)*F(I) 

CONTINUE 

IF CPSO-ACC) 3,3.4 ^ 

PROVIDE PRINTING OF FINAL SOLUTION IF REQUESTED 
IF (IPRINT) 5.6.6 
PRINT 7. MAXC 

FORMATC///5X.'THE FINAL SOLUTION CALCULATED BY NONLIN REQUIRED' 



^ . .. ... . ’F(I)'//(I5.2E17,8)) 

PRINT 9, FSQ 

FORMATC/SX.'THE SUM OF SQUARES IS '.El 7. 8) 

RETURN 

TEST FOR ERROR RETURN BECAUSE FCX) DOES NOT DECREASE 
GO TO CIO. 11, 11. 10. 11). IS 
IF CF5Q*FMIN) 15,20.20 
IF CDD-DSS) 12,12,11 
NTEST = NTEST-l, 

IF CNTEST) 13.14,11 
PRINT 16. NT 

FORMAT C///5X, 'ERROR RETURN FROM NONLIN BECAUSE'. IS, 

1 'CALLS OF CALFUN FAILED TO IMPROVE THE RESIDUALS') 
DO 18 1=1. N 


10 


08200 

#8300 

08400 

08500 

08600 

08700 

08800 

08900 

09000 

09100 

09200 

09300 

09400 

09500 

09600 

09700 

09800 

09900 

10000 

10100 

10200 

10300 

10400 

10500 

10600 

10700 

10800 

10900 

11000 

11100 

il200 

11300 

11400 

11500 

11600 

11700 

11800 

11900 

12000 

12100 

12200 

12300 

12400 

12500 

12600 

12700 

12800 

12900 

13000 

13100 

13200 

13300 

13400 

13500 

13600 

13700 

13800 

13900 

14000 

14100 

T4200 


FROM NONLIN BECAUSE FCX) 
A NEW JACOBIAN*) 


MAXFim GALES 


XCI) = W(NX+I) 

F(I) = WCNF+I) 

CONTINUE 
FSO = FMIN 
GOTO 3 

ERROR RETURN BECAUSE A NEW JACOBIAN IS UNSUCCESSFUL 
PRINT 19 

FORMAT (///5X, "ERROR RETURN FROM NONLIN BECAUSE 
I'FAILED TO DECREASE USING A NEW JACOBIAN*) 

GOTO 17 
NTEST = NT 

TEST WHETHER THERE HAVE BEEN MAXFUN CALLS OF C 
IF CMAXFUN^MAXC) 21,21,22 
PRINT 23,MAXC 

FORMAT(///SX, 'ERROR RETURN FROM NONLIN BECAUSE 
1 15, 'CALLS OF CALFUN*) 

IF (FSO-FMIN) 3,17,17 
PROVIDE PRINTING IF REQUESTED 
IF (IPRINT) 24,24,25 
PRINT 26,MAXC. 

F0RMAT(///5X, *AT THE*,I5, 'TH CALL OF CALFUN WE 
PRINT 8,(I,X(I),FCI),I=1,N) 

PRINT 9,rs6 

GOTO (27, 28, 29/87, 30), IS 

STORE THE RESULT Of THE INITIAL CALL OF CALFUN 
FMIN 5 FSO 
DO 31 1=1, N 
W(NX-H) s XCI) 

WCNF+I) = FCI) 


CALFUN 


NONLIN BECAUSE THERE HAVE BEEN' 


HAVE' 


WCNF+I) = FCI) 

CONTINUE 

CALCULATE A NEW JACOBIAN APPROXIMATION 
IC = 0 
IS 3 3 
IC 3 lC+1 

X(IC) 3 XCIO+DSTEP 
GOTO 1 
K 3 IC 
DO 34 IsljN 

WCK) 3 CFC1)-W(NF+1))/DSTEP 

K 3 K+N 

CONTINUE 

XCIC) 3 WCNX+IC) 

IF CIC-N) 33,35,35 

CALCULATE THE INVERSE OF THE JACOBIAN AND SET THE DIRECTION 
MATRIX 
K 3 0 

DO 36 1=1, N 
DO 37 J31,N 
K 3S K'f 1 

AJINVCI,J) 3 W(K) 

WCND+K) 3 0. 

CONTINUE 

W(NDC+K+I) 3 1 

WCNDC+I) = 1,+FLOATCN-I) 

CONTINUE 

CALL MB01B(AJINV,N,N) 

START ITERATION BY PRIDICTING THE DESCENT i NEWTON MINIMA 
05 * 0 . 

DN 3 0, 

SP 3 0. 

DO 39 1=1, N 


WCK) 


DESCENT i NEWTON MINIMA 
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14300 


xca) = 0. 

a 4400 


F(I) = 0, 

1 14500 

1 14600 
14700 


K ^ I 

00^40 J=1,N 

XCI) = XCi )-«(K)»WCNF+J) 

14800 


FCI) B FCl)-AJiriVCI,J)*W(WF+J) 

14900 


K = K+N 

15000 

40 

CONTINUE 

15100 


DS s DS+XCI)»XCI) 

1 15200 


DN B DN+F(I3*FCI) 

15300 


SP 3 SP+X(I)^F(I) 

15400 

39 

CONTINUE 

15500 

C 

TEST WHETHER A NEARBY STATIONARY POINT IS PRIOICTED 

15600 


IF CFMIN^FMIN-DMM^DS) 41,41,42 

15700 


IF 50 THEN RETURN OR REVISE JACOBIAN 

15800 

42 

GOTO (43, 43, 44), IS 

15900 

44 

PRINT 45 

16000 

45 

F0RMAT(///5X,*ERR0RRETURN FROM NONLIN BECAUSE A NEARBY', 

16100 


1 'STATIONARY POINT OF r(X) IS PREDICTED') 

16200 

43 

GOTO 17 

16300 

16400 

16500 

NTEST = 0 

DO 46 1=1, N 

X(I) 3 WCNX+I) 


46 

16600 

CONTINUE 

1 16700 


GOTO 32 . 

, 16800 


TEST whether TO APPLY THE FULL NEWTON CORRECTION 

' 16900 
17000 

41 

IS s 2 

IF (DN-DD) 47,47,48 

DD s AMAXl (DN,DSS) 

17100 

47 

17200 


DS a 0.25*DN 

1-17300 

1 17400 

49 

TINC * t 

IF (DN-dIs) 49,58,58 

1 17500 

IS 3 4 

1 17600 


GOTO 80 

1 12700 

c ^ 

CALCULATE THE LENGTH OF TTHE STEEPEST DESCENT STEP 

! 17800 

48 

K 3 0 

17900 


DMULT 3 0, 

18000 


DO 51 1=1, N 

18100 


DW s 0 • 

18200 


DO 52 J3l,N 

18300 


K 3 K+1 

18400 


DW 3 DW+WCK)^X(J) 

18500 

52 

CONTINUE 

18600 


DMULT 3 DMULT+DW»DW 

18700 

51 

CONTINUE 

18800 

18900 

19000 

19100 

C 

DMULT 3 DS/DMULT 

DS 3 DS*DMULT»DMULT 

TEST WHETHER TO USE THE STEEPEST DESCENT DIRECTION 

IF (DS-DD) 53,54,54 

TEST WHETHER THE INITIAL VALUEOF DD HAS BEEN SET 

19200 

C 

19300 

54 

IF (DD) 55,55,56 

DD 3 AHAXl (DSS,AMIN1(DM,DS)) 

DS 3 DS/CDMULT*DHULT) 

19400 

19500 

55 

19600 


GOTO 41 

19700 

c 

SET THE MULTIPLIER OF THE STEEPEST DESCENT DIRECTION 

19800 

19900 

20000 

20100 

56 

ANMULT 3 0. 

DMULT 3 DMULT»SQRT(DD/DS) 

GOTO 98 

INTERPOLATE BETWEEN THE STEEPEST DESCENT & THE NEWTON DIRECTION 

c 

^0200 

rO3O0 

53 

SP 3 SP*DMULT 

ANMULT 3 (DD«DS)/( CSP-DS)tSQRT(CSP-DD) ♦♦at (DN-DD )*(DD-DS)>) 


m 
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0400 
0500 
0600 
0700 
20800 
20900 
21000 
1100 
1200 
1300 
21400 
1500 
1600 
1700 
1800 
1900 
2000 
2100 
2200 
2300 
2400 
2500 
2600 
2700 
22800 
22900 
3000 
3100 
23200 
23300 
3400 
3500 
23600 
23700 
23800 
23900 
24000 
24100 
24200 
24300 
24400 
4500 
14600 
[4T00 
14800 
14900 
IS 000 
15100 
1S200 
15300 
15400 
15500 
15600 
15700 
15800 
15900 
26000 
26100 
26200 
16300 
^400 


DMULT = DMULT*(1.-ANMULT) ^ 

CALCULATE THE CHANGE IM X AND ITS ANGLE WITH THE FIRST DIRECTION 
DN = 0. 

SP = 0, 

DO 57 1=1, N 

F(I) = DMULT^XCD+ANMULT^FCI) 

DN = DN+F(I)*F(I) . 

SP = SP+F(I)fW(ND+I) 

CONTINUE 
DS = 0,25*DN 

TEST WHETHER AN EXTRA STEP IS NEEDED FOR INDEPENDENCE 
IF (WCNDC+1)-DTE5T) 58,58,59 
IF (SP»SP-DS) 60,58,58 

TAKE THE EXTRA STEP & UPDATE THE DIRECTION MATRIX 
IS = 2 

DO 61 1=1, N ^ 

X(I) = W(NX+I)+DSTEP*W(ND+I) 

WCN0C+I) = W(NDC+I+1)+1. 

CONTINUE 
W(ND) = 1, 

DO 62 1=1, N 
K = ND+I 
SP » W(K) 

DO 63 J=2,N 
W(K) = WCKtN) 

K = K+N 
CONTINUE 
W(K) = SP 
CONTINUE 
GOTO 1 

EXPRESS THE NEW DIRECTION IN TERMS OF THOSE OF THE DI 
MATRIX, AND UPDATE THE COUNTSIN W(NDC+1) ETC. 

SP * 0 . 

K = ND 
DO 64 1=1, N 
XCI) = DW 
DW = 0. 

DO 65 J=1 ,N 
K » K+1 

DW = DW+F(»3)»W(K) 

CONTINUE 
GOTO (68,66) #IS 
W(NDC+I) = W(NDC+1)+1 
SP = SP+DW^DW 
IF (SP-DS) 64,64,67 
IS = 1 
KK = I 
XCI) = DW 
GOTO 69 
X( I ) = DW 

W(NDC+I) = WCNDC+I+l)+l. 

CONTINUE 

W(ND) =1. . 

REORDER the DIRECTIONS 50 THAT KK IS FIRST 
IF (KK-1) 70,70,71 
KS = NDC-l-KK=<'N 
DO 72 1=1, N 
KS = KS+I 
SP * W(K) 

DO 73 J=2,KK 
W(K) = W(K-N) 


NEW DIRECTION IN TERMS OF THOSE OF 
UPDATE THE COUNTSIN W(NDC+1) ETC. 


THE DIRECTION 


XCI) 

GOTO 

XCI) 


WCNDC+I+l)+l. 


THAT KK 


FIRST 
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26500 
36600 
26700 
26800 
26900 
27000 
27100 
27200 
27300 
27400 
27500 
. 27600 
27700 
27800 
27900 
28000 
28100 
28200 
28300 
28400 
28500 
28600 
28700 
28800 
28900 
29000 
29100 
? 29200 
’29300 
' 29400 

Mitt 

I 29700 
29800 
29900 
30000 
. 30100 
: 30200 
30300 
30400 
30500 
30600 
30700 
30800 
30900 
31000 
31100 
31200 
31300 
31400 
31500 
31600 
31700 
31800 
31900 
32000 
32100 
32200 
32300 

f 2400 
2500 


76 

75 


77 

» ' 

'80 


K s K-N 

73 CONTINUE 
W(K) = SP 

72 CONTINUE 

: GENERATE THE NEW ORTHOGONAL DIRECTION MATRIX 

70 DO 74 1=1, N 

W(NW+I) =1 0, 

74 CONTINUE 

SP a X(1)*XC1) 

K a ND 
DO 75 1=2, N 

DS = SQRT(SP*CSP+X(I)»XCI))) 

DW = SP/DS 
DS a X(I)/DS 
SP = SP+X(I)*X(I) 

DO 76 J=1,N 
K = K+1 

WCNW-^J) = W(NW+J3+X(I-13*W(K) 

WCK) = DW»W(K+N)-DS»W(NW+J) 

CONTINUE 
CONTINUE 
SP = l./SQRTCDN) 

DO 77 1=1, N 
K a K+1 

WCK) = sp*r(i) 

CONTINUE 

CALCULATE THE NEXT VECTOR X,AND PREDICT THE RIGHT HAND SIDES 

V'l ? 

DO 78 1=1, N 

X(I) = WCNX+I)+F(I) 

W(NW+I) = W(NF+T) 

DO 79 J=1,N 

W(NW+I) = W(NW+I)+W(K)»FCJ) 

79 CONTINUE 

FNP = FNP+W(NW+I)f»2 
78 CONTINUE 

C CALL CALFUN USING THE NEW VECTOR OF VARIABLES 

GOTO 1 

C „ UPDATE THE STEP SIZE 

27 DMULT = 0,9*FMIN+0.1*FNP-FSO 

IFCDMULT) 82,81,81 
82 DD = AMAX1(DSS,6.254DD) 

TING a 1, 

IF (FSQ-FMIN) 83,28,28 

C TRY THE TEST TO DECIDE WHETHER TO INCREASE THE STEP LENGTH 

81 SP a 0, 

SS = 0. 

DO 84 1=1, N 

SP = SP+ABS(F(I)*(FCI)-W(NW+I))) 

84 CONTINUE 

PJ a l.+DMULT/(SP+SQRT(SP*SP+DMULT*SS)> 

SP a AMINl(4,,TINC,Pa) 

TINC a PJ/5P 

DD a AMIN1(DM,SP*DD) 

GOTO 83 

C IF FCX) IMPROVES STORE THE NEW VALUE OF X 

ll 
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32600 

§2700 

32800 

32900 

33000 

33100 

33200 

33300 

33400 

33500 

33600 

33700 

33800 

33900 

34000 

34100 

34200 

34300 

34400 

34500 

34600 

34700 

34800 

34900 

35000 

35100 

35200 

35300 

35400 

35500 

f5600 

35700 

35800 

35900 

36000 

36100 

36200 

36300 

36400 

36500 

36600 

36700 

36800 

36900 

37000 

37100 

37200 

37300 

37400 

37500 

37600 

37700 

37800 

37900 


88 

« 

'28 

89 


91 

90 


93 


92 


94 

95 


97 

96 


DO 88 1*1, N 
SP = X(I) 

X(I) s WCNX+1) 

WCNX+I) = SP 
SP * FCI) 

FCI) * WCNF+I) 

W(NF+I) a SP 
W(NW+I) a 
CONTINUE 
IF (IS"U 
CALCULATE _ 

00 89 1=1, H 
XCI) a XtI>*WCNX+I) 

FCI) a F(I)-WCNF+I) 

CONTINUE 

UPDATE THE APPROXIMATIONS TO J AND TO AJINV 
K a 0 

DO 90 lai,N 
W(MW+I) = X(I) 

WCNW-H) a FCI) 

DO 91 Jal,N 

WCMN+I) a WCMW+I)-AJINVCI,J)»FCJ) 


-WCNW+I) 

28,28,50 
THE CHANGES 


IN F AND IN X 


K a K+l 
WCNW+I) a 
CONTINUE 
CONTINUE 
SP a 0. 

SS a 0. 

DO 92 1*1 
a 0. 

93 Jal 


DS 

DO 

DS 


WCNW+I)-WCK)»XCJ) 


-N 


,N 


__ a DS+AJINV(J,I)*XCJ) 

CONTINUE 

SP a SP+DS»FCI) 

SS a SS+XCI)»XCI) 

FCI) a DS 
CONTINUE 
DMULT a 1, 

IF CAflSCSP)-0.1*SS) 94,95,95 
DMULT a 0.8 
PJ a DMULT/SS 

PA a DMULT/CDMULT*SP+C1,-DMULT)*SS) 
K a 0 

DO 96 1=1, N 
SP a PJ»WtNW+I) 

SS a PJ4WCMW+1) 

DO 97 0=1, N 
K a K+1 

WCK) a WCK)+SP*X(J) 

AJINV(I,J) a AJINVCI,J)+SS^FC0) 

CONTINUE 

CONTINUE 

GOTO 38 

END 
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00100 
4)0200 
*00300 
00400 
00500 
00600 
00650 
00660 
00670 
00700 
00800 
; 00900 
■ 01000 
01100 
01200 
01300 
01400 
01500 
01600 
01700 
01750 
01800 
01900 
02000 
02100 
02200 
02300 
02400 
02500 
02600 

J 2700 
2800 
02900 
03000 
03100 
03200 
03300 
03400 
03500 
03600 
03700 
03800 
03900 
04000 
04100 
04200 
04300 
04400 
04500 
04600 
04700 
04800 
04900 
05000 
05100 
05200 
05300 
05400 
05500 
05600 
^5700 


THERMODYNAMIC ANALYSIS OF. COPPER SEGREGATION USING FREE 
ENERGY minimisation TECHNIQUE 
DIMENSION FORTClOO) ,DEliG(100) , 

1PERCC100),AA(100,10) .DFLT C 100 ) , YG(90 ) . XC 100 ) 

1 ,ALPHAC100),XIjAMC100),FSUMC100),YSUMf 100) ^PARPCIOO) 

i,Ya00),AB(l0),R(10,l0),PIC21),B8Cl0) 

UFC100),YC(10).Cfl00) 




CHARACTER *6 SYM(IOO) 

CHARACTER *6 SYMBCIOO) 

NGsNUMBER OF GASEOUS SPECIES 
NC=NUMBER OF CONDENSED SPECIES 

NOsTOTAL NUMBER OF CHEMICAL SPECIES CONCEIVABLE 
MMrNUMBER OF ELEMENTS 
DATA NG,NC,MM,IDEBUG/15,6,6,0/ 

DATA (YCI) ,1=1,7 )/0,00052,, 00596,. 0001 ,.005..0003, .0001 ,,0001/ 
DATA CYCI), 1=8, 11)/, 0327,.5o75,o;6o75,, §05197 
DATA (Y(I),I=12,15.)/4»1E»9/ 

DATA (Y(I),1=16,18)/0.1106, .1316, .733/ 

DATA CY(I),I=l,19)/3’)‘.lE-5/ 

DATA (yCI),I=l,7)/.7E-7,.8E-7,.2E-12,0.0l439,.47E-5,,47E-5, .4E' 
1 ,0,05265/ 

DATA CYCI), 1=8, 11)/0.0155,0. 011961, 0.0024398, 0.0063877/ 

DATA CYCI), 1=12, 15)/. 3E-22,.1E-12,.2E-21,.4E-10/ 

DATA CYC I), 1 = 16, 18)/. 128999,. 096283,. 728671/ 


WRITE (11,3) 
FORMATCIHI ,/lX.30X, 
& SYSTEM%29X,50('- 


'THERMODYNAMIC ANALYSIS OF. CU-FE-C-H-CL-0 
*),//1X,FREE ENERGY AND ATOM COEFFTCIENTSs * ) 


NO=22 

NN=NC+NG 

AP=1. 

AR=1.987 

CRIT=l,0E-5 

AT=298.15 

T=1100, 

DO 31 1=1, NO 

READ (5,8) DELGCD.SYMBCI) 

READ (5,2) (AAC1,J),J=1,MM) 

WRITE (11,602) I, SYMBfI),DELG(I) 

S(AA(I,J),J=1,MM) 

WRITE (11,6) NG,NC,MM 

CPO=0,0645 

FO=0. 73587 

CBN=0. 16444 

HYC=l,-CPa-^CBN 

BB(1)=2,*CP0 

BB(2)=2.*FO 

BB(3)=CBN 

BBC4)=HYC 

BB(5)=HYC 

BBC6)=CP0+3,*FQ 

WRITE (11,5) 

FORMATC/, IX, 'START MOLE FRACTIONS OF CU20 ,FE203 ,C ,HCLCG) Pi ') 
WRITE(11,6 o4) CPO,FO,CBN,HYC 
WRITE (117) 

FORMAT(/,ix, 'ELEMENTAL ABUNDANCE OF CU,C,PE,H,CL,OCSEQUEN )9 ' ) 
WRITE(11,604) (BB(J).J=1,MM) 


^5800 
. 705900 


DO 60 J=1,MM 

1 06000 


ABCJIsO.O 

■ 06100 
' 06200 

610 

DO 610 1=1, NN 

AR(J) = ARCj5+AACI.J)l=y(I) 

06300 


WRITE (11,14) 

06400 

14 

F0RMATC/,1X, 'MASS BALANCE CONSTRAINTS =') 

i 06500 
06600 

41 

WRITE (11,41) . „ , 

FORMATC/IX, 'INITIAL GAS COMPN (SATISFYING MASS BALANCE CONS)=s'> 

06700 

611 

WRITE(11,604) (y(I),l=l,N?0 

06800 

F0RMAT(6B12.5) 

06900 

8 

FORMAT(F10,5,A6) 

07000 

2 

F0RMAT(6F5.2) 

07100 

6 , 

FORMAT (/,16,5X,A6,5X,Fl0.5,5X,6F8.n 

F0RMAT(/,1X,8E15.4) 

07200 

604 

07300 

1032 

FORMAT (IX,12E10.2) 

07400 

9 

DO 1 1=1, NO 

07500 


FORT(I)=((DELG(I)*iOOO.)/(AR*T)) 

07600 

1 

CONTINUE 

07700 


WRITE (11,616) T 

07800 

616 

FORMAT(///,45X, ,TEMPARATURE =' ,Fi0.2) 

07900 


WRITE (11,75) 

08000 

75 

FORMATC45X,24('-')) 

08100 


WRITE (11,2)) 

08200 

21 

F0RMAT(/1X, 'THE IF/RT) VALUES FOLLOW =') 

08300 


WRITE (11,604) (FORT(I),I=l,Nn) 

08400 


MAXNT=50 

08500 


MC=NC 

08600 


MGsNG 

,08700 


XBETAsCRIT 

^8800 


NT=1 

^08900 


DO 33 1=1, MG 

09000 


SyM(I)=SyMB(I) 

09100 

33 

yG(i)=y(i) 

09200 

MN=MG+1 

09300 


MMN=MN 

09400 


NNN=NN 

09500 


IF (MC. EQ. 0) GO TO 300 

09600 


DO 34 I=MN,NN 

09700 


SyM(I)=SYMB(I) 

09800 

34 

YC(i)=y(i) 

09900 

300 

YBARsO. 

10000 


DO 50 1=1, MG 

10100 

50 

YBAR=YBAR+yG(I) 

10200 

10300 

r 

... SET HP AND SOLVE MATRIX 


APRES=DLOG(AP) 

10400 


DO 10 1=1, MG 

10500 


C(I)=FORT(I)+APRES 

10600 


FACsyG(I)/YBAR 

10700 


IF (FAC.LT. l.E-38)FAC=l.E-38 

10800 

10 

F(I)=yG(I)*(C(I)+nLOG(FAC)) 

. 10900 

IF(MC. EQ. 0) GO TO 13 

11060 


J=NG 

11100 


DO 11 1=MN,NN 

11200 


J=J+1 

11300 


C(I)=FORT(J) 

11400 

11 

F(I)=y(I)*C(I) 

11500 

13 

FESUMsO.O 

11600 


DO 43 1=1, NN 

41700 

43 

FESUM=FESUM+F(I) 

T1800 


IF (NT, NE. 1) GO TO 51 
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41900 

^12000 


WRITE (11,61) FESUM 

61 

FORMATC/IX, 'FREE ENERGY SUM OF MIXTRE 

12100 

5 1 

DO 30 J=1,MN 

12200 


DO 30 K=!l,HM 

12300 


ASUMsO. 

12400 


DO 20 1=1, MG 

12500 

20 

SUM=AACI,J)»AA(I,K)*YGCI) 

12600 

ASUMaASUM+SUM 

12700 

12800 

30 

330 

RCJ,K)=ASIIM 

DO 330 1=1, MG 

12900 

ALPHACJ)=ALPHA(J)+AA(X,J)*YGCI) 

13000 


km=mm+mc 

13100 


piM=MM+MC + l 

13200 


MMM=MM+1 

13300 


DO 360 J=MMM,NM 

13400 

360 

ALPHA(J)=0 

13500 


DO 401 1=1, NM 

13600 


DO 401 J=l.l 

13700 

13800 

401 

C0EF(I,J)=ALPHACI) 

MMC=MC+1 

13900 


14000 


IFCMC.EO.O) GO TO 406 

14100 


K=MG 

14200 


DO 402 4=2, MMC 

14300 

14400 


K=K+1 

DO 402 1=1, MM 

14500 

14700 

402 

C0EF(I,J)=AA(K,I) 

MCC=MMC+1. 

14800 

406 

44900 


DO 403 1=1, MM 

*5000 


DO 403 J=MCC,NM 

15100 

403 

K=J-MMC 

15200 

COEF CI,J)=R(I,K) 

15300 


DO 407 I=MMM,NM 

15400 

407 

DO 407 J=2,MMC 

15500 

COEFCT,J)=0. 

15600 


IF CMC. EQ. 0) GO TO 408 

15700 


L=NG 

15800 


DO 404 I=MMM,KM 

15900 


K=0 

16000 


LsL+l 

16100 


DO 404 a^MCC,NM 

16200 


K=K+1 

16300 

404 

C0EFCI,J)=AA(L,K) 

16400 

408 

DO 405 I=MM,NM 

16500 


K=0 

16600 


DO 405 J=MCC,NM 

16700 


K=K+1 

16800 

16900 

405 

COEF(I,J)=ALPHAfK) 

17000 


DO 361 K=l,l 

17100 


DO 361 4=1, MN 

17200 


BSUM(4,K)=0. 

17300 

362 

DO 362 1=1, MG 

17400 

8SUM(4,K)=BSUM(4,K)+AA(I.4)*FCI> 

17500 

17600 

361 

BSUM(J,K)=BSUMCJ,K)+BBC4) 

IF(MC.EO.O) GO TO 409 

DO 363 K=l,l 

17700 


17800 


L=MG 

«47900 

18000 


DO 363 4=MMM,KM 

L=L+1 


* »E18,8) 
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4^100 

718200 

363 

BSUMCJ,K)=CCL) 

18300 

409 

CONTINUE 

18400 


DO 364 K=l,l 

18500 


DO 364 J=NM,NM 

18600 


BSUM(J,K)=0,0 

18700 


DO 364 1=1, MG 

18600 

364 

BSUMCJ,K)=BSUMCJ,K)+F(I) 

18900 


DO 333 J=1,NM 

19000 


K=1 

19100 

333 

FB(J)=B5UM(J,K) 

19200 

CONTINUE 

19300 


CALL SIME0(C0EF,F8,NM,KERR,NM) 

19400 


DO 334 J=1,MN 

19500 


K=1 

19600 

19700 

334 

BSUM(a,K)=FB(J) 

XYBAR=BSUM(1,1) 

19800 


IF(MC.EO.O) GO TO 103 

19900 

C 

NEW X’S FOR CONDENSED SPECIES 

20000 


J=MG 

20100 


DO 101 1=2, MMC 

20200 

20300 


J=J + 1 

20400 

101 

X(J)=BSUM(1 ,1) 

20500 

c 

NEW X*S FOR GASEOUS SPECIES 

20600 

103 

a=o 

20700 


DO 102 I=MCC,NM 

20800 


J=J+l 

20900 

102 

PICJlasBSUMCI,!) 

.21000 


DO 60 1=1, MG 

#1100 

*21200 

60 

FSUMCI)=-FtI)+yG(I3*XYBAR 

DO liO 1=1, HG 

21300 


PSUM=0 

21400 


DO 120 J=1.MM 

21500 

120 

PSUM=PSUM+PI(J)*AA(I,J) 

21600 


YSUMCI)=PSUM*YGfI) 

21700 

110 

X(I)=FSUMa3+YSUMfI) 

21800 

C 

LAMBDA AND DIRECTIONAL DERIVATIVE 

21900 


XLAMBD=1. 

22000 

310 

DO 86 1=1, NN 

22100 


DELT(I)=XCI)-Y(I) 

22200 


IFCDELTCD.GE. 0.) GO TO 86 

22300 


XLAM(I)=YCl)/(-DELTCI)) 

22400 


XLAMBD=DMINJ(XLAMBD,XLAM(I)) 

22500 

86 

CONTINUE 

22600 


DEBAR=0, 

22700 

22800 


DO 87 1=1, MG 1 

22900 

87 

DEBAR=DEBAR+DELT(I) ! 

23000 

93 

DFDLsO.O 

23100 

23200 

96 

DO 88 1=1, MG 

FAC= ( YGC I 5 +XLAMBD#DELT( I ) 3 / (yBAR+XLAMBD*DEBAR) 

23300 


IFCFAC.GT,03 GO TO 82 ; 

23400 


XLAMBD=0.05*XLAMBD 

23500 


GO TO 96 

23600 

82 

DFDL=DFDL+DELT( 13 ♦ (C ( 1 3 +DLOG C C YG C I 3 +XLAMBD»DELT (13 3/ ( YBAR+XL AMB 

23700 

23800 

1«DEBAR333 ; 


23900 

88 

CONTINUE I 

e^4000 

7'24100 


IF(MC.E0.03 GO TO 84 

DO 83 I=HN,NN 
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£4200 

^4300 

83 

DFDL=DFDL+C ( I ) ♦DELT ( I ) 

84 

ir(DFDL,,L.t.O) GO TO 89. . 

24400 

IF(XLAMBD,GT.1.F-38) GO TQ 91 

24500 


KX=KX+1 

24600 


GO TO 191 

24700 

91 

XLAMBD=0,54tXLAMBD 

24800 

GO TO 93 

24900 


IFCIDEBUG.NE.O) WRITE ( 11 . 1021 )XLAMBD 

25000 

1021 

F0RMAT(//,1X,7H,LAMBDA=E15.5//) 

IF (XLAMBO, EQ. 0) GO TO 2001 

25100 

191 

25200 



25300 

C 

NEWSY'S FOR GASEOUS &CONDENSED SPECIES 

25400 


DO 76 1=1, NN 

25500 

76 

Y C I ) =Y C T ) +X1. AMBD+DELT C I ) 

25600 

COMTINUE 

25700 

77 

DO 77 1=1, HG 

YGCI)=Y(i5 

25800 

25900 


IF CMC.EQ. 0) GO TO 351 

26000 


J=flG 

26100 


DO 74 I=MN,NN 

26200 


J=J-H 

26300 

26400 

74 

YC(J)=Y(I) 

26500 

351 

BETA=0. 

26600 

DO 85 1=1, NN 

26700 

85 

8ETA=BETA+DAB5(DE1,T(I)3 . 

26800 


IF (BETA. LT, XBETA) GOTO 140 

26900 


IF (NT. GE. MAXNT) GO TO 600 

27000 


NT=NT+1 

£7100 


GO TO 300 

^7200 

27300 

600 

XBETA=XBETA+CHIT 

27400 


MAXNT=MAXNT+50 

27500 


NTaNT+1 

27600 


GO TO 300 

27700 

140 

PSUM=0. 

27800 

DO 131 1=1, MG 

27900 

28000 

131 

PSUM=PSUM+Y(I) 

28100 

132 

DO 132 1=1, MG 

28200 

PARP(I)=(Y(I3/PSUM) 

28300 

78 

SUMYsO.Q 

28400 


DO 370 1=1, NN 

28500 

370 

SUMY=SUMY+Y(I) 

28600 


DO 340 1=1, NN 

28700 

340 

PERC(I)=(y(T)/SUMYJ- 

28800 


IF CMC, EQ. 0) GO TO 2000 

28900 


J=NG 

29000 


DO 350 I=MN,NN 

29100 

350 

J=J+1 

29200 

PERC(J)=PERC(I) 

29300 

2000 

WRITE (11,16) NT . . 

FORMAT(/1X,28HNO. OF ITERATIONS REQUIRED 
WRITE (11,201) XBETA 
F0RMATC/1X,5HBETA=E12.4/) 

WRITE a 1,22) 

29400 

16 

29500 


29600 

201 

29700 


29800 

22 

F0BMAT(/1X, 'FINAL EOBM COMPN. OF CHEMICAL 

29900 


WRITE (11,32) 

30000 


32 

aoloo 

32 

F0RMAT(/6X,1HI,21X,4HY{I)/) 

30200 


DO 44 1=1, NN 


«T5/) 

SPECIES 



20 


^0300 

^0400 

19 

FORMATC1X,I6,3X,A6,3E20.8,/) 

WRITE (11,19) I,SYMB(I),YCI),PERC(I),PARP(I) 

44 

30500 


WRITE (11,61) FESdM 

30600 


DO 615 J=l,MM 

30700 


ABC.J)=0.0 

30800 

615 

DO 615 1=1, NN 

30900 

AB(J)=ABCJ)+AA(I,J)fy(I) 

31000 


WRITECll ,23) - . 

31100 

23 

F0RMATC/,1X, 'MASS CONSTRAINTS CHECK (WHETHER SATISFIED ?) 

31200 


WRITE (11,604) {AB(J),J=1,MM) 

31300 

31400 

24 

WRITE (11,24) 

FORMAT (/iX, 'LAMBDA VALUES FOLLOWS') 

31500 


WRITE (11,604) (PI(J),J=1.MM) 

31.600 


J=15 

31700 


FORT1=0. 

31800 


DO 4 1=1, MM 

31900 

4 

FORTlsFORTl + (PI ( I) ♦AA( J. I ) ) 

32000 


FORT1=FORT(15)-FORT1 

32100 


PRINT ,F0RT1 . 

32200 


IF (FORTl .GT. 0) GO TO 2001 

32300 


WRITE (11,15) 

32400 

15 

FORMAT (//, 'CU20 IS PRESENT IN THE MIXTURE ') 

32500 

2001 

CONTINUE 

32600 


WRITE (11,25) 

32700 

25 

FORMAT (1 HI) 

32800 

STOP 

32900 


END 

33000 


SUBROUTINE SIME0( A .B ,N . KERR , IDIM) 

33100 



J3200 

C 

THIS SUBROUTINE SOLVES A SYSTEM 

#3300 

c 

OF LINEAR EQUATIONS 

*33400 

c 

BY METHOD OF PIVOTAL CONDENSATION 

33500 

c 

ASNXN MATRIX 

33600 

c 

B: ARRAY CONTAINING N CONSTANTS 

33700 

c 

N: ORDER OF THE SYSTEM - 

33800 

c 

KERR:INDI1CAT0R of subroutine 

33900 

c 

ONE IF SYSTEM SINGULAR 

34000 

c 

ZERO OTHERWISE 

34100 

c 

READ A COLUMN BY COLUMN 

34200 


DIMENSION A(IDIM.IDIM),B(IDIM) 

34300 

- 

DIMENSION L(50).M(50) 

34400 


EPS=lE-6 

34500 


EPSOR=EPS»EPS 

34600 


KERRsO 

34700 

c 

CLEAR OUT PERMUTATION VECTORS 

34800 


DO 3 1=1, N 

34900 


M(I)=0 

35000 

3 

LCI)=0 

35100 

c 

LOOP FOR N PIVOT POINTS 

35200 


DO 14 KP=1,N 

35300 

35400 

35500 

35600 

35700 

35800 

4 

P=0. 

PSOR-0. 

DO 7 1=1, N 

DO 7 J=1,N 

IF {H(I))7,4,7 

IF (L(J))7,5,7 

T=A(I,J) 

35900 

S 

36000 


TSQR=T»T 

36100 


IF (TSQR-PSQR)7,7.6 

-J6200 

36300 

6 

P®T 

PSQRsTSQR 



21 


^6400 


KR=I 

V36500 


KC~U 

36600 

7 

CONTINUE 

36700 


IF CPS0R-EP50P)17.17,8 

36800 



36900 

8 

MCKR)sKC 

37000 


DCKOsl 

37100 

C 

DIVIDE key row by PIVOT 

37200 


DO 10 J=1,N 

37300 


IF CL(J))10,9,10 

37400 

9 

ACKR,J)sa(KR,U)/P 

37500 

10 

CONTINUE 

37600 


0KRsB(KR)/P 

37700 


8CKR)=BKR 

37800 

c 

SUBTRACT MULTIPLE OF KE 

37900 

c 

OTFfER ROWS 

38000 


DO 14 1=1, N 

3B100 


IF (1-KR)H,H,11 

38200 

11 

RsACl,KC) 

38300 


DO 13 Jsl,N 

38400 


IF (L(J))13,12,13 

38500 

12 

A(I,J)=ACI,J)-R»A(KR,J3 

38600 

13 

CONTINUE 

38700 


B(I)=B(I)-RYKBR 

38800 

14 

CONTINUE 

38900 

C 

REORDER RESULTS 

39000 


DO 15 1=1, N 

39100 


IP=M(I) 

39200 

15 

A(IP,1)=B{I) 

J9300 


DO 16 1=1, N 


#9400 

39500 

39600 

39700 

39800 

39900 


16 

17 


ROW FROM 


8(I)=ACI,1) 

RETURN 

ERROR ACTION 
KERRS 1 
RETURN 
END 



